PCT 



WORLD INTELLfiCnJAL PROraRTY ORGANIZA'nON 
bitenutioful Bozeai] 




INTERNAITONAL APPUCATTON PUBUSHED UNDER THE PATEm; CX>OPERATI0N TREATY (PCT), 
(51) Imematkiiml Patent Classification ^ : 



COST 9/DO, 9/40 



Al 



(11) Intermtioiial Publication Number: WO 97/32281 

(43) International Publication Date: 4 Sqttcmbcf 1997 (04JD9.97) 



(21) International Application Number: PCrr/US97/O330O 

(22) International FHing Date: 27 February 1997 (27.02.97) 



(30) Priority Data: 

08/607^88 



27 February 1996 (27.02.96) US 



(60) Parent Application or Grant 
(63) Relsied by Continuation 
US 

Hledoo 



08/607388 (CIP) 
27 Fehniaiy 1996 (27.02,96) 



(71) Applicant (for aU daignated States except US}z INTERVAL 

RESEARCH CORPORATION [USAJSJ; Building C, 1801 
Page Mill Road, Palo Alto. CA 94304 (US). 

(72) Inventors; and 

(75) Invaitors^AppUcants (for US onfy): KOLAROV, Kiasimir 
D. [BGrtJS); 2050 Avy Avenue, Menio Park. CA 94025 
(US). LYNOi Wflliam. C [US/US]; 3331 Tliomas Drive. 
Palo Alto, CA 94303 (US). SCHROlffiR. Pteter [DEOJS]; 
111! BUmchc Street #201. Pasadena. CA 91 106 (US). 
SWELDENS. Wim [BEOJS]; Apartment 2B, 23 Euclid 
Avenue, Summit, NJ 07901 (US). 



(74) Agents: WEAVER, Jeffrey, K. et aL; Hickman Beyer & 
Weaver. P.O. Box 61059. Palo Alto. CA 94306 (US). 



(81) De^gnated States: AL. AM. AT. AU, AZ, BA, BB, EG. BR 
BY, CA. CH. CN. CU. CZ. DE, DK. EE. ES, FI, GB. GE. 
HU, IL. IS. JP. ICE, KG. KP. KR. KZ, LC LK. LR. LS, 
LT. LU. LV. MD. MG, MK. MN. MW. MX. NO. NZ, PL. 
PT, RO. RU. SD. SE. SO, SI. SK, TJ. TM. TR. TT. UA 
UG, US, UZ, VN. ARIPO patent (GH, KE, LS. MW. SD, 
SZ. UG). Eurasian patent (AM, AZ, BY. KG. KZ. MD. RU, 
TJ. TM). European patent (AT. BE. CH, DE. DK, ES. H. 
FR. GB, GR. IE. IT, LU. MC, NL, PT. SE). OAPI patent 
(BF. BJ. CF. CG, CI, CM, OA GN. ML, MR. NE. SN. TO. 
TGI 



Published 

With international search report 

Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
amendments. 



(54) Title: WAVELET BASED DATA COMPRESSION 
(57) Abstract 



A technique for coraprcsxion and 
expansion of a function deftivd upon 
an M-dimensional manifold embedded 
in N-dimenstonal space uses a sec- 
ond generation wavelet tiansfomi and a 
modified zerocree bit-encoding scheme. 
Typically, a function is defined upon 
a two^JimensLooal manifoU embedded 
in tfaree-dimensiooal qMue, such as a 
sphere. A geometric base is chosen 
as a coarse initial model of the man- 
ifold. Second genoaiion waveteis for 
the function are calculated using a tri- 
aiigular subdivision scheme in ocder to 
subdivide the geometric base hi onlcr to 
produce a refined triangular meih. Tte 
wavelet coefiicients are defiivd at the 
vertices of tite trian^es in the triai^ii- 
lar mesh. A tree structure is created in 
which each node of the tree structure 
repietents an associated tiiangle of the 
triangular mesh. Each biangle hi die 
mesh is rccmshrely subdivided into four 
subtriangles and each assodated node in 
the tree structure also has four rfifl^ 



VI 650 



V3 662 



71^622 



T1^624 
V2 664 



n^626 

Tiei2 

V3658 




610 



V36S2 



[2 654 



I^t^^l^a^^J"^ ^ «*triangl« Each wavelet cocfficiem defined at a particular vertex in trtangutar mesh is uniquely assigned to 

^ "^"^^ ^^^^^^^^ nexTsSTlevel of SKTC 
^Z^S^^^^^' ^"V vahK:s of the wavelet coefficients are pr««scd 

rr«?rSl^K • ouipwimg Kte mdicativc of significant nodes and fedr descendants. Sign bits and date bits are also oumut. An 
ex^.«i^jmque.^ bits accord^ ^ 
mveise second genenmon wavelet transf^m is used to symhetize Ac orig^ cocmcienis. An 



FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to identify States 
applicatiocs under the PCT. 



AM 


AmwoU 


AT 


Aaffin 


AU 


AusuiUi 


BB 


BstadM 


BK 


BetfluDi 


BF 


BisktBA PkAO 


BG 


Ba)(tm 


BJ 


Benn 


BR 


Baza 


BY 


Betana 


CA 




CF 


Ctstral Afiicai RepMblic 


CC 


000(0 


CH 


SwitnrtKid 


a 




CM 


Cmucjooo 


CN 




CS 


CocfaMlonku 


CZ 


Czech RepaWic 


DB 


Ocuiiuy 


DK 




EK 




B8 




n 






Rnooe 


GA 


0^ 



to the PCT on die from pages 



GB 


United Kngdora 


GE 




ON 


Gaiaea 


GR 


Gfccoe 


HU 




IB 


tielwd 


IT 


Itdy 


JP 


Japn 


KB 




KG 


KyiiyitiD 


KP 


Detnocntic PM^'i Rqnbbc 




of Kocea 


KR 


Repiditic of Korea 


KZ 




U 




LK 


Sri Ltnka 


LR 


Libsa 


LT 




UJ 




LV 




MC 


Monaco 


MD 


Rfpeblk of Mokkxra 


MG 




ML 


Mdi 


MN 




MR 





parni^lets publishing international 



MW 


Milim 


MX 


Mexico 


NE 




NL 




NO 


Norway 


NZ 


NewZcttand 


PL 


FobBd 


FT 


Poniigal 


RO 


Ronnia 


RU 


jtnniiD ftdwMinn 


SD 


Ssfkn 


SB 


Sweden 


SG 




SI 




SK 


Skmkta 


SN 


Scdcsb) 


sz 




TD 


Chad 


TG 




TJ 




XT 


Tiioidad lad Tofaaso 


UA 


Ubaine 


UG 


Ugai^ 


US 


Uniced Swta of America 


UZ 


Uzbekiiua 


VN 


Via Nam 



wo 97/32281 



PCT/US97/03300 



WAVELET BASED DATA COMPRESSION 



FIELD OF THE INVENTION 

The present invention relates genraally to the compression of data defined on complex 
geometries. More specifically, a compression technique that is adapted to compress data defined 
using a wavelet transformation is described. 

BACKGROUND OF THE INVENTION 

Data compression techniques are useful to represent information as accurately as possible 
with a minimum number of bits in order to minimize the amount of data that must be stored or 
transmitted for use in a computer system. Typically, data to be compressed may represent 
various types of information, such as audio and video information. One of the means primarily 
used for compressing data is the removal of redundant information. At the same time, however, 
the remaining infomiation must accurately represent the data that has been compressed. 

Traditional data compression methods work on data defined on simple geometries. 
Typical geometries used in conventional compression methods include: the real line (e.g., audio), 
a rectangle (e.g., images), or a three-dimensional raster (e.g.. video). However, these 
conventional data compression technologies are not suitable for compression of data defined on 
higher-dimensional or complex geometries such as spheres, general polygonal models, etc. 

There arc numerous examples of complex geometries for which it would be useful to 
develop an appropriate data compression technique. One example is Earth data compression. At 
the moment very large data sets of Earth elevation data exists. A typical example is the ETOP05 
data set. This data set represents topographic data over the entire Earth. The data gives the 
elevation (depth) of the Earth from sea level in meters at a resolution of 5 arc minutes at the 
equator. Even with powerful workstations, manipulating, storing, and transmitting such data sets 
can be quite expensive. Other examples of similar data sets arc satelHte views of the face of the 
Earth and planetary data. Another example of a complex geometry that is useftal for modeling 
mformation is the human head. Another relevant application is the interactive multi-resolution 
viewing of complex data sets. When building an interactive viewer for complex manifolds, it 
would be desirable to only store and display the information needed by the user. 

While various techniques arc available in the prior art for compression of the simple 
geometries, no unified technique currendy exists for the compression of data defined on complex 
and/or higher dimensional geometries. Accordingly, the present invention seeks to provide such 
a technique. 
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. . One definition of a vyayelet Js: ">yavelets arc building blocks that can quicWy derorrclme 

data." The wavelet representation of a data set can be found quickly. More precisely, one can 
switch between the original representation of the data and its wavelet representation in a time 
proportional to the size of the data. The fast decoirelauon power of wavelets is one key to 
5 compression applications. 

In general, wavelets are basis functions that represent a given function at multiple levels of 
detaU. Due to their local support in both space and frequency, they are suited for the sparse 
approximations of functions. Locality in space follows from their compact suppcMt, while locality 
in frequency follows from their smoothness and vanishing moments. In the classic wavelet 
10 setting, i.e.. on the real line, wavelets arc defmcd as the dyadic translates and dUates of one 
particular, fixed function. They arc typically built with the aid of a scaling function. Scaling 
functions and wavelets both satisfy refmemcnt relations (or two scale relations). This means tfiat 
a scaling function or wavelet at a ctttain level of resolution (j) can be written as a linear 
combination of scaling basis functions of the same shape but scaled at one level finer (level j + I). 

1 5 Over the last decade wavelets have become an exceedingly powerful and flexible tool for 

computations and data reduction. They offer theoretical characterizations of smoothness, insights 
into the structure of functions and operators, and practical numerical tools thai lead to faster 
computational algorithms. Examples of their use in computer graphics include surface and 
volume illumination computations, curve and surface modeling and animation, among others. 

20 Given the high computational demand in the quest for speed in computer graphics, the increasing 
exploitation of wavelets comes as no surprise. WhUc these and other appUcations can benefit 
greaUy from wavelet techniques, some of these applications can provide new challenges to the 
underiying wavelet technology. One such challenge is the transformation of data defined on 
complex manifolds. In general, an M-dimensional manifold is a space N such that each point in 

25 N has a neighbortiood homeomorphic to an open ball in R to the M power. For example, a two- 
dimensional (2-D) spherical surface embedded in Euclidean three-dimensional (3-D) space is 
such that each small portion of the sphere appears to be the same as a small portion of Euclidean 
2-D space (i.e.. the plane). This fact caused sonae people to believe that the Earth was flaL 

Classic wavelet consmictions define the wavelet funcUons as dyadic translates and dilates 
30 of one fixed function: "the mother wavelet." See, for example. C.K. Chui. *'An Introduction to 
Wavelets," Academic Press, 1992; A. Cohen, I. Daubechies, and J. Fcauveau. "Bi-orthogonal 
Bases of CompacUy Supported Wavelets," Comm. Pure Appl. Math.. 45:485-560. 1992; and I. 
Daubechies, "Orthonormal bases of compacUy supported wavelets." Conun. Pure Appl. Math,. 
41:909-996. 1988. Because translation and dilation only make sense in a flat Euclidean space. 
35 these classical, or first generation wavelets cannot be used for die transfonnation of data defined 
on more complex manifolds. Therefore, one needs to consider the generalization to second 
generation wavelets. These arc wavelets that arc not necessarily the translates and dilates of one 
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function, but that still enjoy all the powerful properties of first generation wavelets, namely basis 
properties, decorrelation power, and fast algorithms. In general, it can be said that most all of the 
widely used classical wavelets may be represented using second generation wavelets, and it is 
believed that all classical wavelets may be so represented. In other words, second generation 
wavelets are a superset of most all classical wavelets. Thus, the term "second generation 
wavelets" as used herein refers not only to the specific second generation wavelets described, but 
also to those classical wavelets that may be represented using second generation wavelets. 

The basic philosophy behind second generation wavelets is to build wavelets with all 
desirable properties (localization, fast transform) adapted to much more general settings dian the 
real line, e.g., wavelets on manifolds. In order to consider wavelets on a surface, one needs a 
construction of wavelets that are addled to a measure on the surface. In the case of the real line 
(and classical constructions) the measure is dr, the usual translation invariant Ubesgue measure. 
For a sphere, the usual area measure can be denoted by dw, for example. Adaptive constructions 
rely on the realization that translation and dilation arc not fundamental to obtain the wavelets with 
15 the desired properties. The notion that a basis function may be written as a finite linear 
combination of basis functions at a finer, more subdivided level, is mamtained and forms the key 
behind the fast transform. The mam difference with the classical wavelets is that the filter 
coefficients of second generation wavelets are not the same throughout, but can change locally to 
refiect the changing (non-translation invariant) nature of the surface and its measure. 

Classical wavelets and the corresponding filters arc constructed with the aid of the Fourier 
transform. The underlying reason is that translation and dilation become algebraic operations 
after the Fourier transform. In the setting of second generation wavelets, translation and dilation 
can no longer be used, and the Fourier transform cannot be used as a constmction tool. "Ah 
alternative construction is provided by the lifting scheme. 



20 



25 



30 



35 



This new concept in the construction of wavelets, the lifting scheme, has been intit)duced 
recenUy. This concept is described in W. Sweldens, "The Ufting Scheme: A Custom-design 
Construction of Biorthogonal Wavelets". Technical Report 1994:7. Industrial Mathematics 
Initiative, Department of Mathematics. University of South CaroUna, 1994, which is incorporated 
herein by reference in its entirety. This Ufting scheme differs from the classical constructions in 
the sense that it docs not need to use Fourier transforms. Therefore, it is an excellent tool for the 
construction of second generation wavelets. In one implementation, triangular subdivision is 
used to model a three-dimensional manifold, such as a sphere. Such a triangular subdivision of a 
sphere 200 is generally illustrated in Figure 2a. In this figure the surface of sphere 200 has been 
subdivided into triangles in the form of a geodesic sphere. Shown is a uiangle 210 that represents 
one of the many triangles on tfie surface of the sphere. Note that triangle 210 itself forms a part 
of a larger triangle that is represented by triangles 21 1, 210, 212 and 213. It should also be noted 
that triangle 210 itself may be ftiither subdivided according to the same scheme. 
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While modeling techniques thai utilize second generation wavelets to represent data on 

complex manifolds have proven useful, working with such data can be extremely difficult since 
the amount of data required for modeling can be extremely large, and presently, there are no 
suitable mechanisms for compressing such data. Accordingly, the development of a 
5 compression arrangement suitable for compressing such data would be highly desirable. 

SUMMARY OF THE INVENTION 

The present invention embodies a compression technique for compressing and 
uncompressing data representing a function defmed upon an M-dimensional manifold embedded 
in N-dimensional space. With an M-dimensional manifold, a geometric base type is used as a 
10 coarse initial model of the manifold and the geometric base is subdivided into cells. Second 
generation wavelet coefficients for the function are generated using a subdivision scheme in order 
to subdivide the geometric base into cells. The coefficients are defmed upon the cells. 

The invention may also compress and uncompress a function defined upon a two- 
dimensional manifold, where the manifold is embedded in three-dimensional space. In one 

15 aspect of compression, second generation wavelets for the fimction are generated using a 
triangular subdivision scheme in order to subdivide the geometric base in order to produce a 
refined triangular mesh. The wavelet coefficients are defined at the vertices of the triangles in the 
triangular mesh. A tree strucmre is created in which each node of the tree structure represents an 
associated triangle of the triangular mesh. Each triangle in the mesh is recursively subdivided into 

20 four subtriangles and each associated node in the tree structure also has four children, which 
correspond to the four subtriangles. Each wavelet coefficient defmed at a particular vertex in the 
triangular mesh is uniquely assigned to a single one of the triangles at a next higher level of ^ 
subdivision, such that each triangle at the next higher level of subdivision has from zero to three 
assigned wavelet coefficients. 

25 In OTder to compress the wavelet coefficients, the values of the wavelet coefficients are 

processed bit plane by bit plane, by stepping through the nodes of the tree structure. A first 
threshold is set, and a stream of bits is output indicative of significant nodes. A significant node 
is defined as a node diat has either an associated wavelet coefficient diat exceeds the designated 
threshold or a descendant node thai has an associated wavelet coefficient that exceeds the 

30 designated threshold. When a particular significant node is determined, the significant children 
are recursively examined for significant descendants and the output bit stream identifies the 
significant descendants and the significant children nodes. After all nodes have been processed in 
this fashion for this bit plane, a new threshold is set and the next bit plane is processed. The 
technique continues in this fashion until all bit planes for the wavelet coefficients have been 

35 processed. 

-4- 



wo 97/32281 PCT/US97/OJ300 



15 



In a ftirther erabodimeni, when the sign bit of a wavelet coefficient is encountered, it is - 
also output in the bit stream. Also, in a refinement step, for each node already detennined to be 
significant, and not newly determined to be significant, the bit corresponding to the current bit 
plane is output for each wavelet coefficient associated with that node. In another embodiment, the 
5 output bit stream may be entropy encoded using a technique such as Huffman encoding, binary 
encoding, arithmetic encoding, run-length encoding, or LZW encoding. 

In one embodiment, the manifold is a sphere, and the geometric base type may be an 
icosahedron. an octahedron, or a tetrahedron. For each of the geometric base types, the 
calculation of second generation wavelets may use any of a linear, quadratic, cubic or buttcifly 
10 lifting scheme. 

In another embodiment a wavelet coefficient is uniquely associated with a particular 
triangle according to the following technique. As a coefficient is defined at the vertex of a given 
triangle, this vertex will be incident upon the midpoint of an edge thai is shared by two larger 
triangles at the next higher level of subdivision. The coefficient will be uniquely associated with 
one of these two triangles. In one embodiment, all triangles of the mesh are parameterized either 
clockwise or counter-clockwise, and die coefficient is uniquely associated with one of the larger 
triangles according to its parametrization. 

The present invention also embodies a technique to perform uncompression of a function 
defined upon a two-dimensional manifold, where the manifold is embedded in three-dimensional 
space, -nie invention may also be generalized and apply to an M-dimensional manifold 
embedded in N-dimensional space. The ftinction may have been compressed using the above 
described compession technique. The uncompression technique defines a tree structure as 
described above where each node of the tree structure represents a triangle of the triangular mesh. 
From zero to three wavelet coefficients may be associated with each triangle in the mesh 
although the values for these coefficients will not be assigned until the data is input and 
uncompressed. 

The uncompression technique sets a first threshold and inputs a stream of significance bits 
mdKative of significant nodes of the tree stiucture. When a particular significant node is 
determined, the input stream of significance bits further includes bits indicative of the significance 
of children nodes of the particular significant node and the input bit stream identifies the 
significant descendants. In addition, also input are a stream of data bits indicative of values for 
second generation wavelet coefficients, the wavelet coefficients being defined on triangles within 
Ihc refined mesh. Once the wavelet coefficients have been defmed, the technique synthesizes the 
data repn»enting the fiuwtion from the wavelet coefficients in order to reconstnict the function 
defined upon the manifold. A synthesis technique may perform the inverse wavelet transform 
that was used to calculate the wavelet coefficients. 



20 



25 



30 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention, together with further advantages thereof, may best be understood by 
reference to the following description taken in conjunction with the accompanying drawings in 
which: 

5 Figure 1 a is a symbolic representation of an image undergoing a wavelet transform. 

Figure lb is a symbolic representation of the same image in a further iteration of a wavelet 
transform. 

Figure Ic is a symbolic representation of the same image in another iteration of a wavelet 
transform. 

10 Figure Id shows a tree strucmrc of wavelet coefficients used to implement the zerotree 

algorithm. 

Figure 2a shows a spherical object whose surface has been subdivided into triangles 
according to a geodesic subdivision scheme. 

Figure 2b shows a triangle from the surface of the sphere that has undergone a first 
1 5 subdivision as part of a wavelet transform. 

Figure 2c shows a triangle from the surface of the sphere that has undergone a second 
subdivision as part of a wavelet transform. 

Figure 3a is a flowchart showing a process by >yhich data representing a function defined 
upon a complex manifold is compressed according to one embodiment of the present invention. 

20 Figure 3b is a flowchart showing a process by which compressed bits in a data file 

representing a function defined upon a complex manifold arc uncompressed according to one 
embodiment of the present invention. 

Figures 4a, 4b and 4c show a flowchart for a process in which wavelet coefficients 
represented in a tree strucmrc are compressed and output to a file. 

25 Figures 5a, 5b and 5c show a flowchart for a process in which compressed data in a file is 

uncompressed into the original wavelet coefficients. 

Figure 6 shows a triangle from the surface of a sphere that is being subdivided in a 
wavelet transform in accordance with an embodiment of the invention. 

Figures 7a, 7b and 7c illustrate alternative subdivision schemes for subdividing a 
30 geometric shape according to other embodiments of the invention. 
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Figure 8 shows a computer-based system in accordance with an embodiment of the 
present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

The present invention relates to computer readable media that include program 
instnictions for performing various computer-implemented operations. The media and program 
instnictions may be those specially designed and constructed for the purposes of the present 
invention, or they may be of the kind welj known and available to those having skill in the 
computer software arts. Examples of computer readable media include, but are not limited to. 
magnetic media such as hard disks, floppy disks, and magnetic tape; opucaJ media such as CD- 
ROM disks; magneto^piical media such as floptical disks; and hardware devices that are 
specially configured to store and perform program instructions, such as read-only memory 
devices (ROM) and random access memory (RAM). Examples of program instructions include 
both machme code, such as produced by a compiler, and files containing higher level code that 
can be executed by the computer using an interpreter. 

Referring next to Figure 8, a typical computer-based system suitable for implementing the 
present invention is diagrammatically illustrated. Shown is a central processing unit 702 (CPU) 
that is coupled to memory devices including read only memory 704 (ROM) and random access 
meniory 706 (RAM). As is well known in the art, ROM 704 acts to transfer data and instructions 
uni-directionally to the CPU and RAM 706 is used typically to transfer data and instructions in a 
biKlirectional manner. A mass memory device 708 is also coupled bi-directionally to CPU 702 
and provides additional data storage capacity. The mass memory device 708 may be used to store 
programs, data and the like and may take the form of a magnetic or paper tape reader or some ^ 
other well-known device (e,g,. CD-ROM). It will be appreciated that the information retained 
withm the mass memory device 708. may. in appropriate cases, be incorporated in standard 
fashion as part of RAM 106 as virtual memory. CPU 702 is also coupled to one or more 
input/output devices 710 (I/O) that include, but are not limited to, devices such as video monitors, 
track balls, mice, keyboards, microphones, touch-sensitive displays, transducer card readers 
magnetic or paper tape readers, tablets, styluses, voice or handwriting recognizers, or other well- 
known input devices such as, of course, other computers. Finally. CPU 702 optionally may be 
coupled to a computer or telecommunications network using a network connection as shown 
generally at 712. The abovc^lescribed devices and materials will be farmliar to those of skill in the 
computer hardware and software arts. 

As discussed in the background portion of this application, data compression techniques 
typically employ three steps in order to efficiently compress data. These include a transformation 
step, a quantization step and an encoding step. Two of the present inventors previously 
developed the lifting scheme described above for modeling data on complex manifolds m 
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present invention seeks to provide a mechanism for compressing the types of data generated by 
the lifting scheme (and other similar schemes that utilize second generation wavelets). From the 
standpoint of compression, the lifting scheme effectively accomplishes the function of the 
transformation step. 

5 Skipping next to encoding step of the general process for compressing data, a bit encoding 

technique known as the zerotree algorithm has been used in the past to encode classical wavelets. 
This zerotree algorithm method to encode wavelet coefficients was first introduced by J. M. 
Shapiro in "Embedded Image Coding Using Zerotrees of Wavelet Coefficients" 
41(12):34453462. TKRK Trans. Signal Process, 1993. This technique has been extended recently 

10 by A. Said and W. A. Peariman in the manuscript "A New Fast and Efficient Image Codec 
Based on Set Partitioning in Hierarchical Trees", submitted to die IEEE Transactions on Circuits 
and Systems for Video Technology. A further explanation of this technique can also be found in 
"Image Compression Using the Spatial-Orientation Tree", A. Said and W. Peariman, IEEE 
Report 0-7803-1254-6/93, 1993. Each of these publications is incorporated herein by reference. 

15 Said-Peariman have proposed a very fast and efficient method to bit encode wavelet 

coefficients resulting from a conventional wavelet transform. This method and its limitations can 
best be understood by referring to Figures la - Id. Figure la is a symbolic representation of an 
image 100 that has undergone a wavelet transform. Typically, in a wavelet transform an image 
such as the square 102 is subdivided iteratively and at each level a coefficient is associated with a 

20 level of the subdivision. For example. Figure la shows at the coarsest level a square 102 that has 
had a predicted coefficient value A 103 associated with the whole face of the square 102. Further 
iterations of the wavelet transform will further subdivide square 102 into smaller and smaller sub- 
squaics until die pixel level is reached. At each level coefficients will be predicted to correspond 
to the divided sub-squares. 

25 Figure lb shows the next level of iteration of a wavelet transform upon the image 100. In 

this figure, square 102 has been further subdivided into four equally sized sub-squares. Note that 
the coefficient A 103 is still associated with the complete face of the square 102 as indicated by its 
presence in the center of square 102. At this level of iteration, square 102 has been further 
subdivided into squares 104, 106, 108 and 110. Using the wavelet transform and previously 

30 calculated coefficients, new coefficients B 105, C 107, D 109 and E 1 1 1 arc now associated with 
the faces of the newly divided sub-squares 104. 106, 108 and 1 10 respectively. In other words, 
coefficient A 103 is a predicted coefficient value for the face of square 102 while coefficient B 
105 is a predicted value for the square 104. A predicted value might represent the gray scale of a 
particular square, or at the lowest level, a pixel. 

35 Figure Ic represents the next level of iteration upon the image 100 using a wavelet 

transform. Figure Ic demonstrates how the square 102 that had been subdivided in Figure lb is 
now further subdivided. For example, square 104 has now been subdivided into four sub- 
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squares. These sub-squares are numbered 1 12, 1 14. 1 16 and 1 18. Associated with the face of 
each of these smaller sub^squares is a new predicted coefficient, namely F 113. G 1 15, H 117, 
and I 119. For example, if square 112 is at the pixel level then it can be said thai wavelet 
coefficient F 1 1 3 is associated with the whole face of square 1 1 2. Likewise, wavelet coefficient B 
105 is associated with die face of square 104 and wavelet coefficient A 103 is associated with the 
face of square 102. This subdivision process would also be repeated for squares 106, 108 and 
1 10, although not shown. 

Figure Id shows a tree structure 120 tha represents an implementation of the zetotree 
method described by Said-Pteariman. This tree strucnire is also referred to as a quad tree or Q- 
Tree in that the nodes represent squares. The zerotree method first arranges all wavelet 
coefficients in the Q-Tree. At the root of the Q-Tree are the coefficients on the coarsest level. 
(Those coefficients at the top level arc also referred to as die scale coefficients.) Thus, in our 
example, coefficient A 103 represents die coarsest level for the square 102, hence, coefficient A 
103 is at the root of the Q-Tree 120. Each coefficient would then have four "children", namely 
the coefficients of the next finer scale of resolution dial live at die same location. Thus, the 
children of node A will be die coefficients B 105, C 107. D 109 and E 1 1 1. In ottier words, die 
children of node A correspond to die next finer subdivision of square 102. As square 102 was 
subdivided into squares 104, 106, 108 and 1 10, die coefficients diat are associated widi each of 
diese subdivided squares form die chUdren nodes of parent node A. Nodes C, D. E would also 
have children corresponding to a fmer subdivision (as for node B). aldiough diese children are not 
shown. 

The encoding technique of Said-Pearlman now relies on die fact Uiat in smooUi regions of 
images die coefficients decay exponentially. This impUes diat if a certain coefficienris below a t 
direshold tfien its chUdren arc very likely to be below die direshold as well. Thus, in diis likely 
case die whole tree of small coefficients may be discarded. For example, if coefficient B is below 
direshold or close to zero, it is very likely diat its children coefficients F. G, H and I will also be 
below direshold as well. 

The Q-Tree 120 of Figure Id also implicidy encodes die physical location of die 
coefficients. For example, because coefficients B, C, D and E are children of coefficient A, it is 
implicidy known dial coefficient B will occupy die top left comer of square 102, coefficient C will 
occupy die top right comer of square 102. and so on: The zerotree encoding algoridim dius 
enables die coefficients to be compressed and transmitted widiout having to also encode each 
individual cooidinaie for each coefficient. 

However, dierc arc some limitations to dus zerotree algoridim if a geometry odier dian 
square subdivision is used. For example, if a subdivision scheme based upon triangles is used 
and it is desired to associate wavelet coefficients widi die vertices of triangles, certain problems 
arc encountered. As described above, it is often desirable to use a triangular subdivision scheme 
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when modeling complex geometries. These problems will now be discussed in rcfercncc to 
Figures 2a, 2b and 2c. 

Figure 2a shows a sphere 200 that is the subject of data compression using a wavelet 
transform. For example, sphere 200 might represent the Earth, anodier planet or a human head. 
In this figure the surface of sphere 200 has been subdivided into triangles in the form a geodesic 
sphere. Other types of subdivision with or without triangles arc also possible. Shown is a 
triangle 210 that represents one of the many triangles on the surface of the sphere. Note that 
triangle 210 itself forms a pan of a larger triangle that is represented by triangles 211. 210, 212 
and 213. h should also be noted that triangle 210 itself may be further subdivided according to 
the same scheme. 

Figure 2b represents a wavelet transformation process in which large triangle 210 is 
further subdivided into smaller triangles 212, 214. 216 and 218. Because triangle 210 is 
subdivided into four smaller triangles, the Q-Tree representation scheme of Figure Id may be 
^licable. For example, a wavelet coefficient representing the face of triangle 210 would be ai 
the root of the tree and the four wavelet coefficients corresponding to the faces of smaller triangles 
212, 214 . 216 and 218 would be the children of this parent wavelet coefficient This would be a 
straightforward representation of triangle 210 if the coefficients arc associated with the faces of 
the triangles. However, wavelet coefficients are not necessarily always associated with the faces 
of shapes such as triangles or squares. 

In certain subdivision schemes (for example, subdivision using triangles) it may be 
desirable to associate the wavelet coefficients not with the face of a given unit but with its vertices. 
For example, in Figure 2b, triangle 210 has wavelet coefficients defined at jts,, vertices, namely 
wavelet coefficients M 250, N 252, and O 254. Thus, triangle 210 and the wavelet coefficients 
M, N and O represent a first level in a wavelet transform. The next or second level in a wavelet 
transform is reprinted by the subdivision of triangle 210 into smaller triangles 212, 214. 216 
and 218. At this next level, wavelet coefficients M, N and O may be used to predict the second 
level coefficients P 256. Q 258 and R 260. 

Tummg now to Figure 2c, a further subdivision of triangle 210 and a limitaUon of the 
zerotree algorithm is shown. In this figure, triangle 210 has been subdivided again. Top triangle 
212 has now been subdivided into four smaller triangles 220, 222, 224 and 226. Wavelet 
coefficients may now be predicted for the new vertices. As wavelet coefficients M, N and O had 
been used to pt«dict the wavelet coefficients P, Q and R, further prediction occurs sinularly. For 
example, wavelet coefficients M, P and Q have been used to predict wavelet coefficients S 262, T 
264 and U 266. Note how the wavelet coefficients S, T and U reside at the vertices of triangle 
224. 
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It is here that a limitation of the zerotree algorithm becomes apparent. For example, if 
triangle 212 is considered to be a node in the Q-Tree it is unclear which wavelet coefficients 
should be associated with this node. And it is also unclear which wavelet coefficients should be 
associated with the children of triangle 212. For example, it can be seen that triangle 212 may 
5 conceivably be associated with any of the wavelet coefficients M, S, T, P. U or Q. Further 
compounding the problem is the fact that triangle 210 cannot be viewed in isolation. For 
example, referring back to Figure 2a it can be seen that the topmost wavelet coefficient M of 
triangle 210 is also shared with triangles 212, 21 1 and three other triangles. Also, if triangle 220 
is one of the four children in the Q-Tree of triangle 212, it is unclear which wavelet coefficients to 

10 associate with triangle 220. Wavelet coefficients M, S and T are possibilities, but coefficients S 
and T reside ai different levels in the subdivision scheme than coefficient M. That is, wavelet 
coefficient M was produced at a first iteration whereas wavelet coefficient S and T were produced 
at the third iteration. This discrepancy can cause problems in the zerotree algorithm. The 
problem is that given a vertex wavelet coefficient in the triangle tree there is no way to uniquely 

15 assign four children coefficients. This problem is caused by the fact that the wavelet coefficients 
reside at the vertices of the mesh and not at the faces. Note how this problem does not manifest 
itself if wavelet coefficients arc associated with the faces of a subdivision. Referring back to 
Figure Ic, note that square 104 is uniquely associated with wavelet coefficient B and its 
subdivided squares are uniquely associated with coefficients F, G, H and I and so on. Another 

20 problem is that some wavelet coefficients have no children but do have grandchildren. Such a 
configuration prevents application of the original zerotree algorithm since there is no direct link 
betv/een a coefficient and its grandchildren. 

Referring next to Figures 3a and 3b, a method of compressing and uncompressing data in 
accordance with one embodiment of the invention will be described. The compression mediod 

25 of Figure 3a will be described first. The procedure begins at step 301. In this step, the user is 
prompted to input various parameters for use by the procedure. These parameters will be 
described more fiilly below with reference to steps 303 to 315, Typically, these parameters, or 
commands, are input by the user via a command line diat is then parsed by the procedure, as in 
the UNIX operating system. Alternatively, these parameters may be specified by the user via a 

30 graphical user interface (GUI), or the parameters may be input via a script file or the like without 
user interaction. In one embodiment, the invention may be implemented in C++ on a UNIX 
platform.. Graphical visualization may be done, for example, by using the GL library on an SGI 
workstation. 

It is possible for the procedure to operate in an interactive mode. If interactive mode is 
35 selected the user is able to specify parameters to control. In the interactive mode, the user would 
be able to interact with the procedure as subdivision is occurring and the function is being 
compressed. For example, in the interactive mode the user would be able to perform subdivision 
itcratively, viewing the results after each iteration, and would also be able to rotate and view the 
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object from a different angle. Other options available to the user in interactive mode are the ability 
to view the tree structure being formed as a result of the wavelet transformation. If the interactive 
mode is selected^ these options are available to the user as the procedure progresses. The steps of 
the procedure themselves are not changed, only the way in which the user can view the results. 

5 In step 303 the data file is identified. The data file contains the data for a function defined 

upon a manifold. It is this function that will be compressed. The data file also contains a 
representation of the two-dimensional (2-D) manifold itself. The manifold that the function is 
defined on may be a simple rectangle embedded In the plane, as in a normal photogr^h, or it 
may be defined upon a more complex 2-D manifold, such as the spherical surface of the Earth, 
10 which caimot be embedded in the plane. The function may be defined in any of a number of 
ways, including the interpolation of photographic or other sensed data. Techniques for this are 
well known in the ait Intuitively, one may think of the function as being "glued" onto the 
manifold. 

In step 305 the level of subdivision desired is identified. The level of subdivision refers to 
15 how many times a triangle in the geometric base type is subdivided into smaller triangles. For 
example, with reference to Figure 2a, the sphere 200 has been subdivided at the first level. A 
second level of subdivision would mean subdividing each of the triangles shown on the sphere 
200 into four smaller triangles, as shown in Figure 2b. The level of subdivision affects the 
quality of the resulting uncompressed image. For example, a high level of subdivision may 
20 increase the quality of the image, but will take longer to process and requires more memory in a 
computer. 

In step 307 the compression factor is identified. The compression factor is the factor by 
which the data representing a manifold is compressed. For example, a compression factor of 
four (4) would compress the data by four times (4x). In other words, a data file of 4 Megabytes 
25 would be compressed into approximately 1 Megabyte using a compression factor of four (4). A 
high compression factor produces a smaller compressed data file but may result in a lower quality 
uncompressed image. 

In step 309 the number of bit plaries is identified. The concept of a bit plane is known in 
the art. For example, for a graphic image, it is common to represent each value, or pixel, by 8 
30 bits. Thus, there would be 8 bit planes, numbered zero through seven. The i-th bit plane would 
ccmtain the i-th bits for each pixel in the image. In another example, an X-ray image may be 
represented by 1 2 bits for each value, resulting in 12 bit planes. 

In step 311 the geometric base type used for modeling the manifold is identified. It is 
35 desirable to choose a geometric base type that approximates the general shape of the object whose 
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image data is being compressed. For example, in compressing a function defined on the surface 
of the Eartht an icosahedron as the geonietric base type works well to approximate the spherical 
surface of the Earth. For a non-uniform shape, such as an image of a cat, a Dyn sphere works 
well. Dyn spheres are known in the an and are well suited for "blobby" type objects. For a 
5 photograph, a rectangle in the plane as the geometric base type would work best Any geometric 
shape as a base type is possible. It would also be possible to define a new geometric type in 
order to model a unique shape. 



In step 313 the function defined on the manifold is identified. The data representation of 
10 the function will be compressed. For example, the function may represent elevation, color or 
gray scale, or any other image characteristics. For an image of the Earth, the function may 
represent the elevation at a given point on the Earth. In other words, this function may be thought 
of as both a surface, and as a scalar valued function giving height (depth) for each point on a 
sphere. Any arbitrary mathematical function is also possible. 

>5 In step 315 the linear algebraic basis used for the lifting scheme is identified. A basis 

may be defined as a method for generating all elements as a linear combination of the basis 
elements. In the context of second generation wavelets, a basis defines a method for generating 
new coefficients from previously calculated coefficients. This basis will be used by the lifting 
scheme (explained more fiilly below) for a prediction scheme. For example, the basis may be 

20 linear, quadratic, cubic or butterfly (Dyn). It is preferable to use a Dyn basis for compressing data 
representing an image of the Earth. The use by the lifting scheme of a particular basis is 
explained in the reference "Spherical Wavelets: Efficientiy Representing Functions on the 
Sphere", by P. Schroder, W. Swcldens and presented to SIGGRAPH on August 8, 1995 , and is 
hereby incoiporated by reference in its entirety. The use of spherical wavelets for use in the 

25 processing of spherical texture maps is described in the reference "Spherical Wavelets: Texture 
Processing", by P. Schroder, W. Sweldens, Rendering Techniques *95. August 1995, published 
by Springer, and is hereby incorporated by reference in its entirety. 

In step 317 a T-Tree is created. This T-Tree (tree of triangles) contains nodes that will 
represent subdivided triangles of the geomeuic base type. These subdivided triangles are created 

30 in the data analysis step 320 below.. The T-Tree is similar to die Q-Tne of Figure Id tfiat 
represents subdivided squares, except that the T-Trec is a tree representing triangles. As each 
triangle is subdivided into four smaller triangles (as shown in Figure 6, to be explained below), 
each node m the T-Trce (except for leaf nodes) will have four children nodes that represent the 
smaller triangles. As such, each node in the T-Tree represents a triangle in the triangular 

35 subdivision scheme. A wide variety of methods are possible for creating this tree structure. By 
way of example, the triangular subdivision scheme shown in the references "Spherical Wavelets: 
Efficientiy Representing Functions on the Sphere" or "Spherical Wavelets: Texnire Processing" 
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(both referenced above) may be used, to perfonn subdivision, and the T-Tree will represent this 
subdivision. A pointer from each node in the T-Tree may be used to reference a data structure 
that contains information about the represented triangle including its associated wavelet 
coefficients. 

5 In step 319 the error norms for the original data are determined. The original data is the 

image data before transformation or compression. Error norms are used to help evaluate the 
quality of the compression. For example, the maximum and the minimum of the pixel values 
are calculated to evaluate the quality. Other error norms calculated measure the difference 
between the original image and the resulting image. Mathematically, this means establishing a 
10 metric for measuring distances between images. Metrics for measuring these error norms arc 
known in the art. 

In step 320 data analysis is performed to calculate the wavelet coefficients. This is a data 
transformation step in which the original data is transformed into a second generation wavelet 
representation, the data then being represented using second generation wavelet coefficients. This 
15 step performs the subdivision desired upon the manifold and calculates the wavelet coefficients. 
A method for performing this transformation is fully described in the reference "Spherical 
Wavelets: Efficiently Representing Functions on the Sphere", noted above. 

In step 323 the data characteristics for the uncompressed data are calculated These 
characteristics include the error norms and are used for the purposes of peak signal to noise ratio 
20 computation. The calculation for the peak signal to noise ratio may be performed in any manner 
as is known in the an. 

In step 326 the compression function is executed. This function uses a tree based scheme 
that is somewhat analogous to that disclosed in Said-Peariman noted above, and is explained 
more fiilly below with reference to Figures 4a, 4b and 4c. The compression function will 
25 compress the wavelet coefficients representing the fiinction defined upon the 2-D manifold. At 
this point, if a compression file is being written to, the file size before and after 
compression may be compared and displayed to the user. In step 33 1 the error norms for the 
compressed data are calculated. These are the same norms as calculated in step 319, but are 
calculated after compression for comparison purposes. Next, the compression procedure ends. 

30 At this point, compression has been completed. If compression has been performed for a 

purpose such as transmitting the data, storing the data or otherwise manipulating the data, then 
that transmission or storage may be performed at this point and the procedure ends. In some 
cases, however, the user nuy wish to proceed directly to the uncompression procedure of Figure 
3b. For example, for research, testing or other purposes, it may be desirable to uncompress the 

35 data and view the resulting image inunediately. In that case, the user would execute the 
uncompression procedure. Alternatively, if the user begins with or receives compressed data, the 
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user may wish to begin with the uncompression procedure and may not need to perform the 
compression procedure. Accordingly, the uncompression procedure of Figure 3b will now be 
described. 

The uncompression procedure of Figure 3b begins at step 352. In step 354, the user 
identifies which data fiJe contains the compressed data. Of course, the data may also be contained 
in a memory of the computer or may be received in real time over a corranunication line. When 
the data had been compressed, certain information useful to compression is saved in the 
compressed data file. During uncompression, this informauon is first read from the data file in 
order to indicate to the uncompression procedure how the data was compressed. This 
information is read in steps 356, 358, 360, 364 and 366. Of course, if compression is 
standardized, and this information is known a priori, the uncompression procedure would be able 
to make assumptions about this information and would not have to read it from a data file. For 
example, if a decision is made that a tetrahedron will be used for the geometric base type, then the 
uncompression procedure need not read this information but may assume a tetrahedron. 

In step 356 the level of subdivision used for compression is read or received. This level 
corresponds to the level that was idenUfied in step 305 of the compression prtjcedurc. In step 358 
the geometric base type used for compression is read or received This type corresponds to the 
type tfiat was identified in step 31 1 of the compression procedure. In step 360 the basis for the 
lifting scheme used for compression is read or received. This basis corresponds to the basis that 
was identified in step 3 15 of the compression procedure. 

In step 362 the G-Trec is created to serve as a place holder for the wavelet coefficients that 
will be read in from the data file using the uncompression step 368 below. This G-Trce may be 
created as described above in step 404 of Figure 4a, which uses the T-Tree as described in step 
317 of Figure 3a. 

In step 364 the absolute value of the maximum coefficient calculated is read or received. 
This value is used, for example, to convert from fixed point to floating point notation. In step 366 
the number of bit planes used for compression is read or received. This number corresponds to 
the number that was identified in step 309 of the compression procedure. After this information 
has been determined, the uncompression function itself is performed. 

In step 368 the uncompression function is executed. This step is fully explained below 
with reference to Figures 5a, 5b and 5c. The uncompression function wiU take a file of bits that 
had been previously written out using the compression function and uncompress this information 
to geneiatc the set of wavelet coefficients that represent the function defined upon the 2-D 
manifold. Once the coefficients have been generated, a quality analysis may be done. Quality 
analysis involves calculating again the error norms for the newly uncompressed data, and 
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comparing to the previously calculated eriDr norms in. order to detennine the quality of the 
compression. 

In step 370 data synthesis is performed on the wavelet coefficients in order to reconstruct 
as best as possible the data representing the original image. This step is the reverse of step 320 of 
5 Figure 3a, calculating the wavelet coefficients. Data synthesis uses a reverse transform upon the 
second generation wavelet coefficients to produce the image data. A method for performing this 
reverse transformation is also described in the above-identified reference "Spherical Wavelets: 
Efficiently Representing Functions on the Sphere". At this point, the image may be displayed for 
evaluation. 

10 . In step 372 the user is asked if visual evaluation is desired. With an image, it is often 
possible to roughly evaluate the quality of compression by viewing the image before 
compression, and then viewing the same image after it has been compressed and then 
uncompressed. If visual evaluation is not desired, then the uncompression procedure ends. If 
visual evaluation is desired, control moves to step 374, the image is displayed and then the 

1 5 procedure ends . 

Figures 4a, 4b and 4c arc a flowchart illustrating one embodiment of the compression 
function step 326 of Figure 3a. The function begins in step 402. The compression function used 
here is analogous to Algorithm II of Said-Pearlman, although it operates upon a geometry tree (or 
G-Tree) instead of upon a traditional T-Trce as explained more fully below. 

20 In step 404 the G-Tree (geometry tree) is initialized using die T-Tree that was produced by 

the create T-Trcc step 317 of Figure 3a. The G-Tree is related to, but distinct from the T-Tree. 
Different schemes for forming a G-Trec arc possible. One such scheme is described in detail 
below with reference to Figure 6. The G-Tree accommodates triangular meshes of arbitrary two- 
dimensional manifolds, (for example, tht sphere), as contrasted with the purely rectangular 

25 system described in Said-Pcarlman. In the G-Tree, the coefficients live at the vertices of triangles, 
and the G-Trce uniquely associates a vertex with one triangle. In the standard T-Tre«, coefficients 
are defined at the vertices of triangles, but it is not immediately obvious how to associate a 
coefficient with a unique triangle. 

In the standard T-Tree coefficients arc defined at the vertices of triangles because the 
30 analysis step 320 that produces the second generation wavelets creates coefficients that arc defined 
at vertices. Alteraatively. however, it would be possible to generate wavelet coefficients that are 
defined on the faces of triangles, and a corresponding standard T-Tree. In this case, it would not 
be necessary to initialize a G-Tree, because a wavelet coefficient defined on the face of a triangle 
may uniquely be associated with that triangle. 
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In the current embodiment, however, the standard T-Tiee represents coefficients that are 
defined at the vertices of triangles and it is not immediately obvious how to associate a coefficient 
with a unique triangle. Most vertices have six triangles incident upon them. U is necessary to 
associate a coefficient with a particular triangle, even though it may appear that a coefficient may 
be associated with more than one triangle. In other words, the G-Tree will use a tree of 
subdivided triangles to form its backbone, and vertices (and hence, coefficients) will be associated 
widi a particular triangle in the G-Tree. In the G-Tree, then, a specific triangle is picked out with 
which to associate a vertex. There wUI be zero to three vertices mapped to each triangle. TTie 
following provides an explanation of how this mapping is done. 

In the original Said-Pearlman reference a single identification was made between sets of 
coefficients and a single coefficient in die following manner. Each coefficient has three indices, 
the first being the level of subdivision, where the top of the tetrahedn)n, i.e.. the coarsest level is 
numbered zcrt) and finer levels have increasing indices. The second and third indices indicate a 
spatial position within a given level along the X and Y axis (in that order). The value X is a 
smooth coefficient and x* r^resents any of a set of three detail coefficients (since we opmte in 
two dimensions the detail coefficients come in sets of 3. i.e.. * = HL. or LH. or HH). Thus, a 
complete image tetrahedron contains a single ^0,0,0) coefficient at the top and detail coefficients 
at aU finer levels down to some fmest level L-1, i.e.. x*(j.W) with j=0 to L- 1 and k and 1 both = 0 
«o2(j)- I. 

Given this set of coefficients, there is a simple recursive nile by which a given coefficient 
stands as a representative for all coefficients living in the sub-tree underneath it: 

MO AO) enumerates xHUO.O.O), xLH(0.0,0), xHH(O.O.O) 

X*QX.\) enumerates x*(j+l^l). X*(j+U2k+I.21). x*Cj+l.2k.21+l), 

X*a+Uk+1.2I+l) 

where stands for one of HL. LH. or HH, which is arbitrary but fixed. 

This enumeration scheme has to be changed to accommodate surface based compression 
schemes in which coefficients live at die vertices of a uiangular mesh of a 2-D manifold. There is 
no simple nile to use a single coefficient to signify both itself and a whole set of descendent 
coefficients. Thus, the technique of the Said-Pearhnan reference has been modified by 
introducing a way to represent a smglc coefficient as distinct from a way to represent a set of 
coefficients reached from a given triangle in the following manner. 

A coefficient is given as a pain 

(T.E) where T is a triangle pointer and E is an edge designator. 
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The coefficient signified by this is the one sitting at the midpoint of the given edge. The 
special vertices at the top of the hierarchy (4 in the case of the tetrahedron, 8 for the octahedron, 
and 12 for icosahedron) are not the midpoint of any edge. They also happen to be the only sites 
where a smooth coefficient is located They are thus handled separately by simply enumerating 
them From then on each coefficient site, or vertex, is the midpoint of some edge. 

The first difficulty arises in that a given vertex, being the midpoint of a particular 
edge, may have two designators: 

(TJB) and (x, E ) 

where T and x arc the triangles that share the conmion edge E. This can also be explained 
by noting that in working with 2-D manifolds an edge is a one-dimensional (ID) sub-manifold 
and each edge is incident upon either one or two triangles. Because of this a unique way is 
needed to resolve the ambiguity as to which is the proper designator for a given coefficient This 
is done by observing that one of the triangles incident on the edge is parameterized ccHmier- 
clockwise, while the other by necessity is parameterized clockwise. Because the 2-D manifold is 
oriented, the base geometry may be parameterized so that every pair of triangles incident upon a 
common edge have opposite parametrization as to clockwise versus counter-clockwise. With the 
subdivision and parametrization schemes described below, this property is retained at all 
subdivision levels. 

A clockwise triangle is one in which vertices VI, V2 and V3 are enumerated in a 
clockwise direction. The proper designator for a given coefficient is the one for which the triangle 
is labeled counter-clockwise. This rule worics correctly throughout the hierarchy except at the 
edges of the coarsest level, which are 6 for the tcu^edron, 12 for the octahedron, and 30 for the 
icosahedron. For these edges the ambiguity is resolved by choosing as the only legal designator 
amongst the two triangles sharing a given such edge, the one triangle with the lower memory 
address. 

Sets of coefficients are now represented by a triangle alone. It represents all triangles that 
are below it in the hierarchy, i.e., each triangle has 4 constituent children, which themselves have 
possibly 4 constituent children and so forth to level L-1. The set of coefficients represented by 
some triangle T is then simply given by. all coefficients that as a legal designator have: 

(X, E) for some X a triangle desccndent of the root T. 

This representation allows the Said-Pearlman technique to work for data that lives on the 
vertices of some triangle mesh. The triangle, vertex and edge labeling scheme will now be 
described in detail with reference to Figure 6. 
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Figure 6 shows a triangle T 610 that represents one of the triangles on the surface of a 
sphere 200 as shown in Figure 2a, according to an embodiment of the present invention. Triangle 
610 may represent a triangle at the coarsest level, or a triangle at a finer level of subdivision. The 
large triangle T has been subdivided into four medium triangles, each of that may be further 
5 subdivided into smaller triangles. Triangle T has vertices VI 650, V2 654 and V3 652, and has 
been subdivided into medium triangles TO 616 (center), Tl 612 (top), T2 618 (lower right) and 
T3 614 (lower left). This next level of subdivision has produced vertices VI 660, V2 656 and V3 
658, that are shared by new triangles TO, Tl, T2 and T3. 

At the next finer level of subdivision, triangle Tl is subdivided as shown (triangles TO, T2 
10 and T3 would also be subdivided in this fashion, although not shown). Triangle Tl has vertices 
VI 650, V2 656 and V3 658, and has been subdivided into small triangles T1:0 624 (center), 
Tl:l 620 (top), Tl:2 626 (lower right) and Tl:3 622 Cower left). This next level of subdivision 
has produced vertices VI 666, V2 664 and V3 662, which arc shared by new triangles TI:0, 
Tl:l,Tl:2andTI:3. 

15 A triangle is implemented in one embodiment as a record containing four pointers to its 

descendent triangles (diat are nil at the finest level), three pointers to vcitices VI. V2, V3, a bit 
field indicating the parametrization of the triangle (counter-clockwise or clockwise), and a name 
T?, where ? = 0, 1, 2, or 3. The top level triangles arc all counter-clockwise and are named TO. 
Thereafter, child triangles are named after the vertex that they retain from their parent and are 

20 either clockwise or counter-clockwise depending on what the order of the vertices is. For 
example, in Figure 6. the largest triangle is indicated by T 610, and it has child triangles TO, Tl, 
T2 and T3. Triangle TO has a clockwise orientation, while triangles Tl, T2 and T3 all have a 
counter-clockwise orientation. The center triangle does not retain any of its parent vertices. An 
edge of a given triangle is named E? where ? = 1, 2 or 3. The number is chosen to match the 

25 vertex number of the vertex across from the given edge. As a consequence. El is between V2 
and V3, E2 is between V3 and VI, and E3 is between VI and V2. In Figure 6, edges El 670, 
E2 674 and E3 672 are shown for triangle TO. These edges arc shared respecdvely with triangles 
Tl. T2 and T3. Given this, the newly innnoduced vertices in the middle of an edge (due to 
subdivision or instancing of children) are named after their parent edge. El generates VI, etc. 

30 For example. Figure 6 shows vertex V 1 666 as the midpoint of edge E 1 670. 

This naming scheme is self-consistent and unambiguous. Note that triangle TO inherits 
the orientation (counter-clockwise/clockwise) from its parent, while all other triangles have the 
reverse orientation. Accordingly, die coefficient residing at vertex VI 666 may be designated as 
(Tl, El) or as (TO, El). As Tl has a counter-clockwise orientation, die corrca designation is 
35 (Tl, El). This designation is arbitrary, and it is contemplated dial the correct designation may be 
the one for which the triangle has a clockwise orientation. Thus, using die above representation. 
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the T-Tree may be converted into a G-Tree in. which .wavelet coefficients residing at vertices arc 
associated with a unique triangle. 

It is possible to utilize other schemes for forming the G-Tree. By way of example, for the 
subdivision scheme of Figure 6, a G-Tree may be formed where the nodes of the G-Tree are any 
5 of triangles, edges, or vertices. The base geometry forms the root of the G-Tree and immediately 
below it are each of the vertices, edges and triangles of the base geometry. To subdivide a triangle 
in this fashion, its edges are first subdivided. Each edge subdivides into two edges and a (central) 
vertex. Each of these is a G-Tree child of the subdivided edge. The triangle is now subdivided 
into four triangles and three edges (the interior edges) that become the G-Tree children of the 
10 subdivided triangle. Each vertex is now uniquely assigned to either the base geometry or to an 
edge. Each edge is uniquely assigned to one of the base geometry, another edge, or a triangle. 
Each triangle is uniquely assigned to one of the base geometry or another triangle. 

The preceding scheme for subdivision may be generalized to M-dimensional manifolds. 
Vertices, edges, and triangles are all instances of cells (0-cells. 1 -cells and 2-cells respectively). 

15 Cells have a dimension, and cells of dimension M are called M-cells. This concept of cells is well 
known to those of skill in the ait A generalization to M-dimensional manifolds may be 
performed as follows. Subdividing a cell requires a previous subdivision of the cells on its 
boundary (all of which are of strictly lower dimension). The subdivision of a 0-cell is itself. The 
cells in the subdivision become the G-Tree children of the subdivided cell, and are associated with 

20 the subdivided cell. The cells of the base geometry are the G-Tree children of the base geometry 
itself, and each is associated with the base geontetiy. 

^ Within this generalized scheme there are a number of methods for performing the actual 
cell subdivision. One such subdivision scheme is barycemric subdivision which is known in the 
art Figure 7a illustrates the barycentric subdivision of a triangle. Baiycentric subdivision 

25 subdivides an M-cell into a (central) 0-cell and a K+1 ceU for each K-cell of the subdivided 
boundary of the M-cell. Each of the K-cells in the subdivision have K+1 0-cells on their 
boundary and is a K-dimensional generalization of the triangle (called a simplex^ a 2-simplex is a 
triangle, a 3-simplex is a pyramid). Note that the barycentric subdivision of a triangle has 13 cells 
(the central 0-cell. six 1 -cells from the six vertices of the subdivided boundary, and six 2-cells 

30 from the six edges of the subdivided boundary). 

Another scheme for the actual ceU subdivision is bi-conic subdivision. Bt-conic 
subdivision of a triangle is shown in Figure 7b, and of a quadrilateral in Figure 7c. Bi-conic 
subdivision subdivides an M-cell into a (central) 0-cell and an M-K cell for each K-cell of the un- 
subdivided bcmndary of the M-cell. Each of the K-cells in the subdivision have at least K+1 0- 
35 cells on their boundary. A 2-ccll in the subdivision will be a quadrilateral. These various 
subdivision schemes are illustrative of the ways in which a cell may be subdivided. Once the G- 
Tree has been created here in step 404. the compression fimcdon continues. 
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In step 406 the wavelet coefficients arc converted to binary form. Specifically, each 
wavelet coefficient in scientific floating point notation is converted to fixed point binary signed 
magnitude notation. This conversion may be accomplished by walking the G-Tree and 
converting each coefficient in turn. 

In step 408 the significance bits "S(N)" for each wavelet coefficient are calculated In 
general, the calculation of significance bits is known in the ait This significance function is 
analogous to the one described in Said-Pearlman except that the tree on which the S(N) bits are 
calculated is not die Q-Tree of Said-Pearlman but the G-Tree. The "S" function is a mask that 
indicates whether a bit plane is significant for a sub-tree. More specifically, for a given wavelet 
coefficient of N bits, the corresponding significance function S(N) will be a if a bit at bit 
position N is not a leading *'0'' (or if S(N) for a descendant coefficient is a "1"), and S(N) will be 
a **0'* if the bit at position N is a leading zero and S(N) for every descendent coefficient is also a 
**0*'. For example, if a wavelet coefficient has a binary representation of " 0001 0100 then the 
corresponding S(N) function would be "0001 1 1 11 the first three binary digits of S(N) being 
each a "0*^ because the wavelet coefficient has three leading zeros. 

As alluded to above, S(N) for a given wavelet coefficient in the G-Tree will also be a " 1 " 
if for that bit position N, there exists a descendant wavelet coefficient with a corresponding S(N) 
that also has a "1" in that position. In other words, S(N) for a given wavelet coefficient **W" in 
the G-Tree is the logical "or** of the S(N) functions corresponding to all wavelet coefficients that 
20 are descendants of **W** in the G-Tree. Using the above example, if **W** is a wavelet coefficient 
represented as " 0001 0100 and there exists a descendent of "X" with a coefficient 
represented as " 001 1 0000 then S(N) for "X" is " OOU 1111 and S(N) for "W* would 
also be "0011- 1111"; 

Steps 410 through 466 of Figures 4a, 4b and 4c implement a procedure analogous to 
25 Algorithm II of Said-Pearlman, although the data being processed is quite different since a G- 
Tree is used This procedure uses various lists and sets to facilitate its implementation. LIS 
stands for "List of Insignificant Sets", LIN is "List of Insignificant Nodes" and LSN is "List of 
Significam Nodes.** The sets used are 0(T), all offspring or children of a node T, D(T), all 
descendants of a node T. and L(T) which is D(T) minus 0(T) (all descendants of node T except 
30 for immediate children). This procedure will now be explained. 

In step 410 the bit position N is set equal to the most significant bit of the binary 
representation of the wavelet coefficient For example, if 8 bits are used to represent the wavelet 
coefficients, then N = 7. For the least significant bit, N = 0. This bit position N is also preferably 
output to a file, that will contain the compressed data when the compression step 326 of Figure 3a 
35 is finished The other "output" steps referred to below also indicate that the bit is to be output to 
this same file. However, it is not necessary that the bits be output to a file. It is also 
contemplated that the compressed bits may be stored in memory, written to a disk, or otherwise 

-21- 



10 



wo 97/32281 



PCr/US97/0J300 



stored in a computer system. Furthermore, the bits need not be stored locally, but may be 
transmitted upon output via a computer network, telephone line, microwave, satellite transmission 
or the like to be processed at a remote location. 

In step 41 2 the ordered lists LIS (Ust of Insignificant Sets) and UN (List of Insignificant 
5 Nodes) are initialized to contain all of the top level nodes. A top level node represents a wavelet 
coefficient at the coarsest level of the subdivision, as a result of the analysis step 320 of Figure 3a. 
In step 4 14 the list LSN (List of Significant Nodes) is initialized to the empty UsL 

Step 416 represents a looping structure in which index T is used to step through each 
entry in the UN. At first, T is initialized to the first entry in the UN. If the end of the LIN has 
10 been reached then control transfers to step 426 of Figure 4b. If the end of the UN has not been 
reached then the loop continues through steps 418 to 424. 

In Step 418 the significance bit, S(N), is output for node T. Because there may be from 0 
to 3 wavelet coefficients associated with this node T, there will be from 0 to 3 significance bits 
ouqjut in this step, depending upon how many coefficients arc associated with this node. Step 

15 420 tests whether the significance of the bit at bit position N for any of the coefficients associated 
with node T is equal to I. If not, control rcmms to step 416 and T is set equal to the next entry in 
the UN. If so, step 422 is executed in which node T is moved to the LSN. Next, in step 424, the 
sign of each of the coefficients at node T is output After step 424 control returns to step 416 in 
which T is set equal to the next entry in the LIN. As soon as the end of the LIN is reached control 

20 moves to step 426 of Figure 4b. 

Step 426 of Figure 4b is a looping stmcture that steps through each entry in the US until 
the end of the LIS is reached. Initially, T is set equal to the first entry in the LIS. When the end 
of the LIS is reached control moves to step 458 of Figure 4c. Step 428 tests whether node T is 
marked as representing set D. If so, then steps 430-448 are executed. If not, then T is marked as 

25 representing set L and control moves to step 450 of Figure 4c where steps 450-456 are executed. 
In step 430 all of the significance bits at bit position N for all the coefficients in set LXJ) arc 
(HiQnit Step 432 tests whether the significance bit at bit position N for any of the coefficients in 
set D(T) is equal to one. If the significance bit is not one then control returns to step 426 and the 
counter T is set equal to the next entry in the US. If the significance bit is a one then control 

30 nwves to step 434. 

Step 434 is a looping structure in which an index U is used to step through each entry in 
the set 0(T). Initially U is set equal to the first entry in set 0(T). If the end of set 0(D is reached 
then control moves to step 446. In step 436 the significance bit at bit position N for each of the 
coefficients at node U is output Step 438 tests whether the significance bit at bit position N for 
35 any of the coefficients associated with node U is equal to one. If not, control moves to step 444. 
If so, control continues to step 440. In step 440 node U is added to the LSN. In step 442 the sign 
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of each of the coefficients at node U is output. After step 442 control returns to step 434 and U is 
set equal to the next entry in set 0(T). In step 444 node U is added to the LIN and then control 
returns to step 434. Step 446 tests whether set UJ) is empty. If so. then control returns to step 
426 and T is set equal to the next entry in the LIS. If set L(T) is not empty then in step 448 node 
5 T is moved to the end of the LIS and is marked as representing set L. After step 448 control 
remms to step 426. 

In step 450 of Figure 4c the significance bits at bit position N for each of the coefficients 
in the set L(T) are output Step 452 tests whether the significance bit at bit position N for any of 
the coefficients associated with the set L(T) is equal to a one. If the significance bit is not a one 
10 then control returns to step 426 of Figure 4b. If the significance bit is a one then step 454 is 
executed. In step 454 each member of set 0(T) is added to the end of the LIS and is marked as 
representing set D. In step 456 node T is removed from the LIS and control then returns to step 
426 of Figure 4b. 

Step 458 represents a looping structure in which an index T is used to step through each 
15 entry in the LSN. Initially T is set equal to the first entry in the LSN. If th& end of the LSN is 
reached then control moves to step 464. Step 460 tests whedier the entry T in LSN has been 
included in the most recent sorting pass thixiugh LIN and LIS. If so then control returns to step 
458 and T is set equal to the next entry in LSN. If not, then step 462 is executed. In step 462 the 
Nth most significant bit of the absolute value of each of the coefficients at node T is output. After 
20 step 462 control returns to step 458. In step 464 the counter N is decremented by one. Step 466 
checks whether N is equal to zero. As N represents the number of bit planes for the coefficients, 
if N is equal to zero then step 326 of Figure 3b is done. If N is not zero then control returns to 
■ step 4 16 of Figure 4a and T is set equal to the next entry in the UN. 

Turning now to Figures 5a, 5b and 5c, step 368 of Figure 3b "execute uncompression 
25 function" will now be explained in detail. Steps 510 through 566 of Figures 5a, 5b and 5c 
implement an uncompression procedure analogous to Algorithm II of Said-Peariman, although 
the following procedure uncompresses data that has previously been compressed, whereas a 
literal reading of Said-Pearlman produces a compression algorithm. Nonetheless, the reader will 
appreciate that this uncompression procedure is analogous. Figure 5a shows a flowchart 
30 lepresenling the uncompression procedure step 368 of Figure 3b. The procedure begins at step 
502. 

In step 510 the most significant bit N in the binary representation of the wavelet 
coefficients is read in finm the output file. N now represents the current bit position, starting with 
the most significant bit This bit position N is preferably input from a file that contains die 
35 compressed data after die compression step 326 of Figure 3a has finished. The other "input" 
steps referred to below also indicate that the bit is to be input from this same file. However, it is 
not necessary that the bits be input from a file. As described above, it is also contemplated dial 
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the compressed bits may be read from memory, a disk, or other storage, area in a computer 
system. Furthermore, the bits need not be read locally, but may be read from a computer 
network, telephone line, microwave, satellite transmission or the like and then be processed 
locally. 

5 In step 512 the ordered lists LIS and LIN are initialized to contain all of the top level 

nodes. These top level nodes may be read from the ou^t file where they were written 
previously in the compression step. A top level node represents a scale coefficient at the coarsest 
level of the subdivision, as a result of the analysis step 320 of Figure 3a. In step 5 14 the list LSN 
is initialized to the empty list 

10 Step 516 represents a looping structure in which index T is used to step through each 

entry in the UN. At first, T is initialized to the first entry in the UN. If the end of the UN has 
been reached then control transfers to step 526 of Figure 5b. If the end of the UN has not been 
reached then the loop continues through steps 5 1 8 to 524. 

In Step 518 the significance bit, S(N), is input for each of the coefficients associated with 
15 node T. Step 520 tests whether the significance of the bit at bit position N for any of the 
coefficients of node T is equal to I. If not, control returns to step 516 and T is set equal to the 
next entry in the UN. If so, step 522 is executed m which node T is moved to the LSN. Next, in 
step 524, the sign of each of the coefficients at node T is input After step 524 control returns to 
step 516 in which T is set equal to the next entry in the UN. As soon as the end of die UN is 
20 reached control moves to step 526 of Figure 5b. 

Step 526 of Figure 5b is a looping structure that steps through each entry in the US until 
the end of the LIS is reached. Initially, T is set equal to the first entry in the LIS. When the end 
of the US is reached control moves to step 558 of Figure 5c. Step 528 tests whether node T is 
marked as representing set D. If so, then steps 530-548 are executed. If not, then T is marked as 

25 representing set L and control moves to step 550 of Figure 5c where steps 550-556 arc executed. 
In step 530 all of the significance bits at bit position N for all the coefficients in set D(T) are input 
Step 532 tests whether the significance bit at bit position N for any of the coefficients in set D(T) 
is equal to one. If the significance bit is not one then control returns to step 526 and the counter T 
is set equal to the next entry in the US. If the significance bit is a one then control moves to step 

30 534. 

Step 534 is a looping stiucmrc in which an index U is used to step through each entry in 
the set 0(T). Initially U is set equal to the first entry in set 0(T). If the end of set 0(T) is reached 
then control moves to step 546. In step 536 the significance bit at bit position N for each of the 
coefficients at node U is input Step 538 tests whether the significance bit al bit position N for 
35 any of the coefficients associated with node U is equal to one. If not, control moves to step 544. 
If so, control continues to step 540. In step 540 node U is added to the LSN. In step 542 the sign 
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of each of the coefficients at node U is input After step 542 controi returns to step 534 and U is. 
set equal to the next entry in set 6(V). in step 544 node U is added to the LIN and then control 
returns to step 534. Step 546 tests whether set UJ) is empty. If so, then control rctums to step 
526 and T is set equal to the next entry in the US. If set L(D is not empty then in step 548 node 
T is moved to the end of the LIS and is marked as representing set L. After step 548 control 
returns to step 526. 

In step 550 of Figure 5c the significance bits at bit position N for the coefficients in the set 
UD are input Step 552 tests whether the significance bit at bit position N for any of the 
coefficients in the set UJ) is equal to a one. If die significance bit is not a one then control returns 
to step 526 of Figure 5b. If the significance bit is a one then step 554 is executed. In step 554 
each member of set 0(T) is added to the end of the US and is mariced as representing set D. In 
step 556 node T is removed from die US and control dien rctums to step 526 of Figure 5b, 

Step 558 represents a looping strucnirc in which an index T is used to step through each 
entry in the LSN. Initially T is set equal to the first entry in the LSN. If die end of the LSN is 
reached dien control moves to step 564. Step 560 tests whether the entry T in LSN has been 
included in the most recent sorting pass tfirough LIN and US. If so then control renims to step 
558 and T is set equal to die next entry in LSN. If not dien step 562 is executed. In step 562 die 
Ndi most significant bit of die absohite value of each of die coefficients at node T is input After 
step 562 control returns to step 558. In step 564 die counter N is decremented by one. Step 566 
checks whedier N is equal to zero. As N represents die number of bit planes for die coefficients, 
if N is equal to zero dien step 368 of Figure 3b is done. If N is not zero dien control remms tci 
step 516 of Figure 5a and T is set equal to die next entry in die LIN. 

Aldiough die foregoing invention has been described in some detail for purposes of clarity 
of understanding, it will be apparent dial certain changes and modifications may be practiced 
widiin die scope of die appended claims. For instance, die image or manifold on which a 
function is to be compressed need not be a sphere, but may be any odicr two^mensional 
manifold. Also, die geometric base type may be any geometric shape, or even a specially defined 
geometry. The invention may also be practiced upon an M-dimensional manifold diat is 
embedded in N-dimensional space, where M and N are integers and N is greater dian M. In 
addition, die subdivision scheme may use any geometric sh^ for performing die subdivision 
and need not be restricted to triangles. By way of example, in higher dimensions die geometric 
base may be subdivided into cells. A cell is a unit of subdivision and is known in die art For a 
2-D manifold dje cell may be any polygon, and for a 3-D manifold die cell may be any polytope. 
And although a second generation wavelet transform is described to perform die data analysis and 
generate die wavelet coefficients, it is contemplated diat odier data transformation techniques may 
be used to generate die coefficients and build up a tree struchire. 
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And within the G-Tree, although only one technique has been described for uniquely 
associating a coefficient defined at a vertex with a triangle, it may be possible to use other 
techniques for associating a coefficient with one unique triangle. Additionally, a coefficient may 
be associated with the face of a triangle, meaning that the coefficient would be uniquely associated 
5 with that triangle. Also, although an embodiment has been described that operates upon a static 
data set defined upon a sphere, for example, the invention may also be used to compress and 
uncompress dynamically-changing, tinie-dependent data sets. Such data sets occur in 
applications such as animation, modeling of namral phenomena, etc. 

EXEMPLARY EMBODIMENT 

10 In one embodiment of the invention, optimization of the present invention for high- 

dimensional compression is achieved through combining discrete second generation wavelet 
transforms with an extension of the embedded zerotree coding method. Other optimizations 
include the use of arithmetic coding, different scaling and norms of the wavelet coefficients. 
Also, different prediction schemes in the zerotree mechanism may be used. Combining these 

1 5 techniques can lead to a significant improvement in compression performance. A discussion of 
other embodiments arc described in Xompression of Functions Defined on Surfaces of 3D 
Objects", K. Kolarov and W. Lynch. Proceedings of the Data Compression Conference (1997), 
which is hereby incorporated by reference for jurisdictions permitting such incorporation. 

This embodiment discusses optimizations through the use of (1) Scaling. (2) Parameter 
20 Variation, and (3) Entropy Coding which are discussed below. 

Example data sets used for optimization are an Earth image data set and the Lena image. 
For the Earth example a function that is the elevation (with respect to sea level) of the Earth is 
compressed. This topographic function is approximated from the ETOPOlO data set which 
samples the Earth from satellites every 10 arc minutes (1.5 million points approximately 1 1 miles 
25 apart). For the base complex, an icosahedron is used and is subdivided eight times, yielding 

655362 vertices and 1310720 triangles. The vertices of this complex are radially projected onto a 
concentric geoid (vertices average about a 22 arc minute separation - about 25 miles). The 
topographic elevation at each vertex is then determined by interpolation of the ETOPOlO data set 
(2:1 reduction in points). 

30 For comparison with existing data compression techniques the Lena image is compressed 

using an embodiment of the present invention (the function being compressed is the gray-scale 
level at each pixel). A base complex includes two triangles, thus forming a square. This complex 
is subdivided nine times yielding a complex with 5 1 3^ vertices and 2 • 5 1 2^ triangles. The 
function value at each of the 513^ vertices is deteimined by interpolation from a 5 12^ gray-scale 

35 Lena image. 



-26- 



wo 97/32281 



PCT/US97/03300 



The evaluation of a compression technique may be perfonned in a wide variety of ways, 
and may depend upon the compression technique used to return a meaningful result By way of 
example, for the present invention, evaluation of compression by looking at the amount of 
infomialion that is read during decompression in relation to the amount of information that gets 
5 written out by the compression module worics well. Parameters that may be analyzed to 

determine compression performance include bits per pixel, a PSNR, and a compression ratio. 

One of the input parameters to an embodiment of the present invention is the number of 
bit planes used for representing each wavelet coefficient. Variation of die number of bit planes can 
affect compression performance. During compression the wavelet coefficients generated by the 

1 0 lifting are scaled and then ftirther processed by the zerotree construction . The resulting output 
stream consists of coefficients (in order of magnitude) from the most significant to the least 
significant bit planes. This output stream inchides all of the bit planes for the coefficients as 
truncated by the zerotree hypothesis. At decompression it can be decided (based on time and 
space resources) to read in all or just a portion of the output bit stream written during 

1 5 compression. This ability reflects the "progressive coding" character of the code. 

A value for bits per pixel is a good measure of compression performance which together 
with the PSNR can be used for comparison with other algoriUims. The notion of bits per pixel 
corresponds to bits per (wavelet) coefficient, since the coefficients arc the ones that represent the 
data structure. Using the Lena image as an example, for 10 bit planes read in during 
20 decompression, the actual size of the file is 369 1 89 bits for the Lena image. Those bits 
correspond to 513^= 263169 wavelet coefficients, yielding a rate of 1.4029 bits per pixel. 

A compression ratio is a widely used characteristic for comparison of algorithms. The 
compression ratio is calculated with respect to the amount of data needed to write out directly all 
of the wavelet coefficients. Each coefficient takes 8 bits for complete representation. Thus, in the 
25 example of the Lena image, 263 169 coefficients wiU need 2105352 bits for complete 

representation. Using 10 bit planes per coefficient the present invention reads (and first writes) 
369189 bits as outlined. Thus, the true compression ratio is 6: 1 for these parameters. A 1 : 1 ratio 
can achieved by taking more bit planes per coefficient, for example, 15 bit planes achieves that 
ratio. 

30 SCALING 

During decompression the wavelets arc descaled by the inverse of the scaling factor used 
during compression. As explained above, optimization for the different norms can be achieved 
by scaling the coefficients after the transform analysis and before the zerotree compression for a 
given threshold. Multiplication by unity optimizes for the norm. To optimize for the 
35 norm the wavelet coefficients may be scaled by 2"^ where j is the subdivision level at which the 
coefficient is attached This basically corresponds to shifting the bits in the coefficients' 
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representation one to the left To optimize for the norm, the wavelet coefficients may be scaled 

by 4"'. This corresponds to two bits shift in the coefficients. For the general norm the scaling 

zi 

factor is 4 . The contrast sensitivity curve for the human visual system is best approximated by 
using the norm and its associated scaling. 

5 Comparative results for the Lena image with and L_ reveal that the performance 

curves for the L, and the scaling are actually the same. All points using the I, norm lie on the 
curve for but are much earlier on the curve than the corresponding bit planes points from 
norm. L, points further along the curve can be obtained by using more bit planes for the 
coefficients (e.g. 15 bit planes yield a point with 32.36 dB for 0.45 b/vertex). Since 1, vs. 

10 involves shifting two positions left for vs. one position left for L^^ it is natural to expea that a 
snrudler output stream results, and yields a smaller bits per pixel rate and a coirespondingly 
smaller PSNR for the same number of bit planes read. Also, any scaling of L, or performs 
better than no scaling at all (i.e. norm). Using scaling results in a 6 dB improvement over the 
whole range of the performance curve. scaling of coefficients is used for the below variations 

1 5 in parameters and coding for this embodiment 

PARAMETER VARIATION 

Varying the different parameters used in the technique of the present invention affects 
compression performance. This section considers the influence of the number of levels of 
subdivision used, the role of the number of bit planes used to represent the coefficients, and the 
20 use of different prediction schemes. 

The levels of subdivision of the base manifold affects compression performance. There is 
an advantage to using as many subdivision levels as possible in order to cover all of the data 
points in the original database. Adding one more level of subdivision increases the PSNR with 
more than 5 dB for both the Lena image and the Earth data. Thus, as much subdivision as 
25 possible is advantageous, although it may be constrained by hardware limitations. 

The number of bit planes to represent the wavelet coefficients also affects compression 
performance. A wide variety of number of bit planes may be used. Ten bit planes works well, 
although in some cases it may be desirable to have more bit planes in order to obtain a higher 
PSNR and better visual perception quality. Naturally, for images like Lena the maximal bits per 
30 pixel (coefficient) rale that is needed is 8 bits per pixel, since that is what it takes to represent all 
coefficients explicitly. To obtain a 1 : 1 basis for visual comparison of images, it is useful to 
determine the maximal number of bit planes needed. 

For example, the Lena image using 9 levels of subdivision can obtain a 1 : 1 compression 
ratio with a 15 bit plane representation of the coefficients. The resulting PSNR is 66.56 dB for 
35 8. 15 bits per pixel. Calculating the performance curve for reading in at decompression using 15, 
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13 2, 1 bit planes of the 15 bit planes written out at compression yields the same curve 

with the new points beyond 10 bit planes lying further to the right on the curve. 

Using a linear vs. a butterfly prediction scheme may also affects compression 
performance. The above examples use the linear prediction scheme as outlined above. In the 
linear prediction scheme, the new wavelet coefFicients generated at the midpoints of any given 
edge arc calculated as the mean of the two coefficients associated with the vertices of the edge. 
Another scheme is the butterfly scheme where the new coefficient is a linear combination of eight 
neighboring vertices from the previous level, i.e. a larger stencil is used for prediction of new 
coefficients. Using the Earth data example, 8 levels of subdivision and butterfly prediction brings 
about a 4% improvement over the linear case. The geometry of the mesh approximation of the 
Earth sphere being quite regular niay be responsible for this slight increase. Improvements in 
performance may be seen when using larger stencils for high-dimensional objects with sharp 
edges and concavities. 

ENTROPY CODING 

The addition of an entropy coder and modifications to the coding scheme also affect the 
performance of the present invention. 

Arithmetic coding of the output binary stream can be advantageous. The result from 
applying a wavelet transfrara, scaling of the wavelet coefficients and the zetotree compression is 
an output bit stream. The performance curve can be further enhanced by applying Huffman 
coding or arithmetic coding to this output stream. 

A wide variety of arithmetic coders may be used. By way of example, the wcH-known 
Witten coder may be used. In one embodiment, the present invention generics the arithmetic 
coded output stream during decompression in paraUel with reading in the compressed data. The 
PSNR is not affected by the arithmetic coding; the size of the bit stream is smaller. Using the 
Lena image as an example, using 9 levels of subdivision reveals that the arithmetic coder brings 
on the average 33% in^vement over straight binary coding. 

An alternative way of buOding the G-tree in the zerotree implementation helps to achieve a 
significant improvement of more than 12% compared to the results without this alternative way. 
One implementation of the zcnrtrec compression using G-trecs is described in detail above. The 
main structure of that implementation bears simUarity with the SPIHT (Algorithm 0) described 
in "A New Fast and Efficient Image Codec Based on Set Partitioning in Hierarchical Trees^' 
referenced above. This implementation uses new data stmctures (the G-tree) and control flow to 
account for the more complicated topology of the spaces to be processed. Varying the 
initialization and the branching phases affects compression performance as wUl now be described. 
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G rr is defined as a list of triangles in the G-tree, GIV and GS V are lists of vertioES in the 
G-trce. The first part of the algorithm for initialization is: 

Set n = N = sig(D(R)) output N. 

Initialize GIT to the triangles in the base complex (as type A triangles). 
Initialize GIV to the vertices in the base complex. 
Initialize CSV to empty. 
Here D(R} is the set of descendants of the root of the tree R. 

In other embodiments, the top level coefficients at the vertices (4 for the square defming 
the Lena image and 20 for the icosahedron defming the Earth approximation) may be written out 
diiBctly in the output stream. This creates a choice for the initialization of the data structures. 

A first option is to initialize GIT with the top level triangular faces (2 for the square and 20 
for the icosahedron) and initialize GIV to empty (since the lop level vertices are not processed). 
Using this option the fan out on the next level per vertex is maximum 4 and good depth can be 
achieved in the tree. The significance of a vertex is the value that is calculated for iL The 
significance of a triangle is based on the significance of the triangles into which it is subdivided as 
well as the significance of its vertices. 

A second option is to initialize GIV with the midpoints of the edges of the base triangles 
in err (5 in the case of the square and 30 for the icosahedron). The original top level vertices are 
not used because they art not>occssed. This strucmrcs calculates significance ofall triangles of 
Type A by looking at their children (triangles) in the tree and their vertices. The significance of all 
triangles of type B (already encountered) is calculated using the significance of the descendants 
except itself and its children. The fan out of this approach is 8 per node and the G tree is generally 
not as deep as in the previous case. 

The first option has very similar performance curve to the second option. The 
improvement for the second option is about 4% on the average for the Lena image and 1 .75% on 
the average for the Earth. However if arithmetic coding is added to the wider G-tree generation a 
significant improvement is achieved. The compression curve becomes on the average 10% better 
than the resolution. The improvement in performance is more than I dB. 

Thus, the butterfly scheme allows a 4% improvement over the linear one. For a span of 4 
chUdicn per vertex the arithmetic coding improves the perfonnance by 4% on the average. For the 
span of 8 children the arithmetic coder brings about 10% improvement Good results are 
achieved when as many subdivision levels as necessary are used to cover the original data. More 
bit planes allow compression ratios closer to I : I and help to obtain base picmres for comparison. 
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Also, it is beneficial to use scaling and works well as an appropriate norm for a good 
compression curve. Using a combination of these techniques brings about 12% improvement in 
performance with respect to the results without these techniques. 

The invention has appUcaUons in telecommunications, on the Internet, in database and 
computer graphic applications and in similar industries. Any application where complex data is 
being manipulated, compressed or transmitted may benefit from use of the invenuon. Therefore, 
the described embodiment should be taken as Ulustrative and not restricUve, and the invention 
should be limited solely by the scope of the following claims. 
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CLAIMS 

1. A computer-implemented method for compressing data representing a function 
defined upon an M-dimensional manifold, the manifold being embedded in N-dimensional space 
where M and N arc integere and N is greater than M, a geometric base being used as a coarse 
5 initial model of the manifold, the method comprising the steps of: 

calculating second generation wavelets for the function using a subdivision scheme for 
subdividing the geometric base into cells in order to produce a refined mesh, the wavelet 
coefficients being defined on cells within the refined mesh; 

defining a tree structure in which each node of the tree structure represents an associated 
10 cell of the refined mesh, there being a plurality of levels in the tree, the tree levels including a lop 
level tha includes at least one top level node, wherein each tree level below the top level includes a 
plurality of nodes that are children of parent nodes in a next higher level, wherein chUdrcn nodes 
associated with a particular parent node represent ceUs of the subdivided cell associated with the 
particular parent node, wherein each wavelet coefficient at a particular level in the refined mesh is 
1 5 uniquely defmed on a particular cell in the rcfmed mesh at that level and is uniquely assigned to a 
single one of the ceUs. such that each cell may have assigned to it any number of wavelet 
coefficients; 

setting a first threshold and outputting a stream of bits indicative of significant nodes, 
wherein a significant node is defmed as a node that has at least one of, associated wavelet 

20 coefficient that exceeds die designated threshold and descendant node that has an associated 
wavelet coefficient tfiat exceeds the designated threshold, wherein when a particular significant 
node is determined, the output stream of bits hirther includes bits indicative of die significance of 
children nodes of the particular significant node and die significant chUdrcn are recursively 
examined for significant descendants and the output bit stream identifies the significant 

25 descendants. 

2. A method as recited in claim 1 wherein each tunc a significant node is identified, 
an indicia of the sign of each significant coefficient associated with the significant node is also 
output 

3. A method as recited in claim 1 further comprising a refinement step that refmes 
30 the setting and outputung step, die refinement step including the substcps of: 

setting a current direshold to a next threshold; 
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outputting a stream of bits indicative of newly significant nodes at the current 
threshold, wherein a newly significant node is defined as a node that is significant at the current 
threshold but was not significant at any previous threshold; and 

each time a new coefficient that exceeds the current threshold is identified, outputting an 
5 indicia of the sign of the new coefficient 

4. A method as recited in claim 3 further comprising the step of repealing the 
refinement step wherein the current threshold is set to a threshold that is one-half of (he previous 
threshold. 

5. A method as recited in any of claims 1 -4 wherein the M-dimensional manifold 
10 embedded in N-dimensional space is a 2-dimensional manifold embedded in 3-dimensional 

space. 

6. A method as recited in any of claims 1-5 wherein the subdivision scheme is a 
triangular subdivision scheme, the cells of the subdivided geometric base are triangles and the 
refined mesh is a refined triangular mesh that contemplates that each base triangle in a selected 

1 5 level of the triangular mesh is divided into four subtriangles at the next tower level of subdivision, 
such that the vertices of the subtriangles include the vertices of the base triangle and the midpoints 
of the edges of the base triangle. 

7. A method as recited in any of claims 1-6 wherein the calculating step uses a lifting 
scheme with a basis that is selected fh}m the group consisting of linear, quadratic, cubic and 

20 butterfly. 

8. A method as recited in claim 6 wherein each wavelet coefficient is defined at the 
vertex of a triangle within the refined triangular mesh and the particular vertex at which each 
wavelet coefficient is defined is incident upon the midpoint of an edge belonging to two larger 
triangles at a next higher level of subdivision, and the wavelet coefficient is uniquely associated 

25 with one of tfte two larger triangles. 

9. A computer program product comprising a computer-usable medium having 
computer-readable code embodied thereon for compressing data representing a function defined 
upon an M-dimensional manifold, the manifold being embedded in N-dimensional space where 
M and N are integers and N is greater than M, a geometric base being used as a coarse initial 

30 model of the manifold, the computer program product comprising the following computer- 
readable program code for effecting actions in the computer 

program code for calculating second generation wavelets for the function using a 
subdivision scheme for subdividing the geometric base into cells in order to produce a refined 
mesh, the wavelet coefficients being defined on cells within the refined mesh; 
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program code for defining a tree structure in which each node of the tree structure 
represents an associated cell of the refined mesh, there being a plurality of levels in the tree, the 
tree levels including a top level that includes at least one top level node, whoem each tree level 
below the top level includes a plurality of nodes that are children of parent nodes in a next higher 
5 level, wherein children nodes associated with a particular parent node represent cells of the 
subdivided cell associated with the particular parent node, wherein each wavelet coefficient at a 
particular level in the refined niesh is uniquely defined on a particular cell in the rcfmed mesh at 
that level and is uniquely assigned to a single one of the cells, such that each cell may have 
assigned to it any number of wavelet coefficients; 

10 program code for setting a first threshold and outputting a stream of bits indicative of 

significant nodes, wherein a significant node is defined as a node that has at least one of, 
associated wavelet coefficient that exceeds the designated threshold and descendant node that has 
an associated wavelet coefficient that exceeds the designated threshold, wherein when a particular 
significant node is determined, the output stream of bits further includes bits indicative of the 

1 5 significance of children nodes of the particular significant node and the significant children are 
recursively examined for significant descendants and the output bit stream identifies the 
significant descendants. 

10. A computer-implemented method for compressing data representing a function 
defined upon a two-dimensional manifold, the manifold being embedded in three-dimensional 
20 space, a geometric base being used as a coarse initial model of the manifold, the method 
comprising the steps of: 

calculating second generation wavelets for the function using a triangular subdivision 
scheme for subdividing the geometric base in order to produce a refmed triangular mesh, the 
wavelet coefficients being defined at the vertices of the triangles in the triangular mesh; 

25 defining a tree structure in which each node of the tree structure represents an associated 

triangle of the uiangular mesh, there being a plurality of levels in the tree, the tree levels including 
a top level that includes at least one top level node, wherein each tree level below the top level 
includes a plurality of nodes that are children of parent nodes in a next higher level, wherein 
children nodes associated with a particular parent node represent triangles of the subdivided 

30 triangle associated with the particular parent node, wherein each wavelet coefficient at. a particular 
level in the triangular mesh is uniquely defined at a particular vertex in the triangular mesh at that 
level and is uniquely assigned to a single one of the triangles at a next higher level of subdivision, 
such that each triangle at the next higher level of subdivision has in the range of zero to three 
assigned wavelet coefficients; 

35 setting a first threshold and outputting a stream of bits indicative of significant nodes, 

wherein a significant node is defined as a node that has at least one of, associated wavelet 
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coefficient that exceeds the designated threshold and descendant node that has an associated 
wavelet coefficient that exceeds the designated threshold, wherein when a particular significant 
node is determined, the output stream of bits further includes bits indicative of the significance of 
children nodes of the particular significant node and the significant children are recursively 
examined for significant descendants and the output bit stream identifies the significant 
descendants. 

11. A method as recited in claim 10 wherein each time a significant node is identified, 
an indicia of the sign of each significant coefficient associated with the significant node is also 
output 

12. A method as recited in claim 10 further comprising a refinement step that refines 
the setting and outputling step, die refinement step including die subsieps of: 

setting a current threshold to a next threshold; 

Quitting a stream of bits indicative of newly significant nodes at the current 
threshold, wherein a newly significant node is defined as a node that is significant at the current 
threshold but was not significant at any previous threshold; and 

each time a new coefficient diat exceeds the current direshold is identified, outputting an 
indicia of the sign of the new coefficient. 

13. A method as recited in claim 12 further comprising the step of repeating the 
refinement step wherein the current direshold is set to a direshold that is one-half of die previous 
direshold. 

14. A mediod as recited in any of claims 10-13 wherein die manifold is a sphere and 
die geometric base is selected from die group consisting of an icosahedron, an octahedron and a 
tetrahedron. 

15. A mediod as recited in any of claims 10-14 wherein die triangular subdivision 
scheme contemplates tiiat each base triangle in a selected level of die mesh is divided into four 
subtriangles at die next lower level of subdivision, such diat die vertices of die subtriangles 
include die vertices of die base triangle and die midpoints of die edges of die base triangle. 

16. A mediod as recited in any of claims 10-15 wherein die calculating step uses a 
lifting scheme widi a basis dim is selected from die group consisting of Unear, quadratic, cubic 
and butterfly. 

17. A mediod as recited in any of claims 10-16 wherein die particular vertex at which 
each wavelet coefficient is defined is incident upon die midpoint of an edge belonging to two 
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larger triangles at a next higher level of subdivision, and the wavelet coefficieni is uniquely 
associated with one of the two larger triangles. 

1 8. A computer program product comprising a computer-usable medium having 
computer-readable code embodied thereon for compressing data representing a function defined 
upon a two-dimensional manifold, the manifold being embedded in three-dimensional space, a 
geometric base being used as a coarse initial model of the manifold, the computer program 
product comprising the following computer-readable program code for effecting actions in the 
computer 

program code for calculating second generation wavelets for the function using a 
triangular subdivision scheme for subdividing the geomeuic base in order to produce a refmed 
triangular mesh, the wavelet coefficients being defined at die vertices of the triangles in the 
triangular mesh; 

program code for defining a tree structure in which each node of the tree structure 
represents an associated triangle of the triangular mesh, there being a plurality of levels in the tree, 
the tree levels including a top level that includes at least one top level node, wherein each tree level 
below the top level includes a plurality of nodes that are children of parent nodes in a next higher 
level, wherein children nodes associated with a particular parent node represent triangles of the 
subdivided triangle associated with the particular parent node, wherein each wavelet coefficient at 
a particular level in the triangular mesh is uniquely defined at a particular vertex in the triangular 
mesh at that level and is uniquely assigned to a single one of the triangles at a next higher level of 
subdivision, such that each triangle at the next higher level of subdivision has in the range of zero 
to three assigtjed wavelet coefficients; . 

program code for setting a first threshold and outputting a stream of bits indicative of 
significant nodes, wherein a significant node is defined as a node that has at least one of, 
associated wavelet coefficient diat exceeds the designated threshold and descendant node that has 
an associated wavelet coefficient that exceeds the designated threshold, wherein when a particular 
significant node is determined, the output stream of bits further includes bits indicative of the 
significance of children nodes of the particular significant node and the significant children are 
recursively examined for significant descendants and the output bit stream identifies the 
significant descendants. 

19. A computer-implemented method for uncompressing data representing a function 
defined upon an M-dimensional manifold, the manifold being embedded in N-dimensional space 
where M and N arc integers and N is greater than M, a geometric base being used as a coarse 
initial model of the manifold, the geometric base being subdivided into cells using a subdivision 
scheme to produce a refined mesh, the method comprising the steps of: 
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defining a tree stmcture in which eadi node of the tree structure represents an associated 
ceU of the refined mesh, each node airanged to accommodate associated second generation 
wavelet coefficients, there being a plurality of levels in the tree, the tree levels including a top Jevel 
that includes at least one top level node, wherein each tree level below the top level includes a 
plurality of nodes that arc children of parent nodes in a next higher level, wherein children nodes 
associated with a particular parent node represent cells of a subdivided cell associated with the 
particular parent node; 

setting a first threshold and inputting a stream of significance bits indicative of significant 
nodes, wherein a significant node is defined as a node that has at least one of, associated wavelet 
coefficient that exceeds the designated threshold and descendant node that has an associated 
wavelet coefficient that exceeds the designated threshold, wherein when a particular significant 
node is determined, the input stream of significance bits further includes bits indicative of the 
significance of children nodes of the particular significant node and the input bit stream identifies 
the significant descendants; 

inputting a stream of data bits indicative of values for second generation wavelet 
coefficients, the wavelet coefficients being defined on cells within the refined mesh, wherein each 
wavelet coefficient at a particular level in the refined mesh is uniquely defined on a particular cell 
in the refined mesh at that Jevel and is uniquely assigned to a single one of the cells, such that each 
cell may have assigned to it any number of wavelet coefficients; and 

synthesizing the data representing die function from die wavelet coefficients in order to 
reconstruct the function defined upon the M-dimensional manifold. 

20. A mediod as recited in claim 19 wherein each time a significant node is identified, 
an indicia of the sign of each significant coefficient associated with the significant node is also 
input 

21. A method as recited in claim 19 hiither comprising a refinement step that refines 
die setting and inputting steps, die refinement step including die substeps of: 

setting a current direshold to a next dircshold; 

inputting a stream of bits indicative of newly significant nodes at die current 
threshold, wherein a newly significant node is defined as a node diat is significant at die current 
direshold but was not significant at any previous threshold; 

inputting a stream of data bits corresponding to die current direshold, die data bits 
indicative of values for die second generation wavelet coefficients. 
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22. A method as recited in claim 21 further comprising the step of repeating the 
refmement step wherein the current threshold is set to a threshold that is one-half of the previous 
threshold. 

23. A method as recited in any of claims 19-22 wherein the M-dimensionaJ manifold 
5 embedded in N-dimensional space is a 2-dimensional manifold embedded in 3-dimensional 

space. 

24. A method as recited in any of claims 19-23 wherein the subdivision scheme is a 
triangular subdivision scheme, the cells of the subdivided geometric base are triangles and the 
refined mesh is a refined triangular mesh that contemplates that each base triangle in a selected 

1 0 level of the triangular mesh is divided into four subtriangles at the next lower level of subdivision, 
such that the vertices of the subtriangles include the vertices of the base triangle and the midpoints 
of the edges of the base triangle. 

25 . A method as recited in any of claims 1 9-24 wherein the synthesizing step uses a 
lifting scheme with a basis that is selected from the group consisting of linear, quadratic, cubic 

15 and butterfly. 

26. A method as recited in claim 24 wherein each wavelet coefficient is defined at the 
vertex of a triangle within the refined triangular mesh and the particular vertex at which each 
wavelet coefficient is defined is incident upon the midpoint of an edge belonging to two larger 
triangles at a next higher level of subdivision, and the wavelet coefficient is uniquely associated 

20 with one of the two larger triangles. 

27. A computer program product comprising a computer-usable medium having 
computer-readable code embodied thereon for uncompressing data representing a function 
defined upon an M-dimcnsional manifold, the manifold being embedded in N-dimensional space 
where M and N arc integers and N is greater than M. a geometric base being used as a coarse 

25 initial model of the manifold, the geometric base being subdivided into cells using a subdivision 
scheme to produce a refined mesh, the computer program product comprising the following 
computer-readable program code for effecting actions in the computer: 

program code for defining a tree structtire in which each node of the tree structure 
represents an associated cell of the refined mesh, each node arranged to accommodate associated 
30 second generation waveirt coefficients, there being a plurality of levels in the tree, the tree levels 
including a top level that includes at least one top level node, wherein each tree level below the top 
level includes a phirality of nodes that are chUdren of parent nodes in a next higher level, wherein 
children nodes associated with a particular parent node represent celU of a subdivided ceU 
associated with the particular parent node: 
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program code for setting a first threshold and inputting a stream of significance bits 
indicative of significant nodes, wherein a significant node is defined as a node that has at least one 
of, associated wavelet coefficient that exceeds the designated threshold and descendant node that 
has an associated wavelet coefficient that exceeds the designated threshold, wherein when a 
5 particular significant node is determined, the input stream of significance bits ftirther includes bits 
indicative of the significance of children nodes of the particular significant node and the input bit 
stream identifies the significant descendants; 

program code for inputting a stream of data bits indicative of values for second generation 
wavelet coefficients, the wavelet coefficients being defined on cells within the refined mesh, 
10 wherein each wavelet coefficient at a particular level in the refmed mesh is uniquely defmed on a 
particular cell in the refmed mesh at that level and is uniquely assigned to a single one of the cells, 
such that each cell may have assigned to it any number of wavelet coefficients; and 

program code for synthesizing the data representing the function from the wavelet 
coefficients in order to reconstruct the function defined upon the M-dimensional manifold. 

15 28. A computer-implemented method for uncompressing data representing a function 

defined upon a two-dimensional manifold, the manifold being embedded in three-dimensional 
space, a geometric base being used as a coarse initial nxxiel of the manifold, the geometric base 
being subdivided using a triangular subdivision scheme to produce a refined triangular mesh, the 
method comprising the steps of: 

20 defining a tree structure in which each node of the tree structure represents an associated 

triangle of the triangular mesh, each node airanged to accommodate associated second generation 
wavelet coefficients, there being a plurality of levels in the tree, the tree levels including a top level 
that includes at least one tc^ level node, wherein each tree level below the top level includes a 
plurality of nodes that arc children of parent nodes in a next higher level, wherein children nodes 

25 associated with a particular parent node represent triangles of a subdivided triangle associated with 
the particular parent node; 

setting a first threshold and inputting a stream of significance bits indicative of significant 
nodes, wherein a significant node is defined as a node that has a least one of, associated wavelet 
coefficient that exceeds the designated threshold and descendant node that has an associated 
30 wavelet coefficient that exceeds the designated threshold, wherein when a particular significant 
node is determined, the input stream of significance bits further includes bits indicative of the 
significance of children nodes of the particular significant node and the input bit stream identifies 
the significant descendants; 

inputting a stream of data bits indicative of values for second generation wavelet 
35 coefficients, the wavelet coefficients being defined at the vertices of the triangles in the refmed 
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triangular mesh, wherein each wavelet coefficient at a particular level in the triangular mesh is 
uniquely defined at a particular vertex in the triangular mesh at that level and is uniquely assigned 
to a single one of the triangles at a next higher level of subdivision, such that each triangle at the 
next higher level of subdivision has in the range of zero to three assigned wavelet coefTicients; and 

5 synthesizing the data representing the function from the wavelet coefficients in order to 

reconstruct the function defined upon the two-dimensional manifold. 

29. A method as recited in claim 28 wherein each time a significant node is identified, 
an indicia of the sign of each significant coefficient associated with the significant node is also 
input. 

10 30. A method as recited in claim 28 further comprising a reHnement step that refines 

the setting and inputting steps, the refinement step including the substeps of: 

setting a current threshold to a next threshold; 

inputting a stream of bits indicative of newly significant nodes at the current 
threshold, wherein a newly significant node is defined as a node that is significant at the current 
1 5 threshold but was not significant at any previous threshold; 

inputting a stream of data bits corresponding to the current threshold, the data bits 
indicative of values for the second generation wavelet coefficients. 

31. A method as recited in claim 30 further comprising the step of repeating the 
refinement step wherein the current threshold is set to a threshold that is one-half of the previous 

20 threshold. 

32. A method as recited in any of claims 28-31 wherein the manifold is a sphere and 
the geometric base is selected from the group consisting of an icosahedron, an octahedron and a 
tetrahedron. 

33. A method as recited in any of claims 28-32 wherein the triangular subdivision 
25 scheme contemplates that each base triangle in a selected level of the triangular mesh is divided 

into four subtriangles at the next lower level of subdivision, such that the vertices of the 
subtriangles include the vertices of the base triangle and the midpoints of the edges of the base 
triangle. 

34. A method as recited in any of claims 28-33 wherein the synthesizing step uses a 
30 lifting scheme with a basis that is selected from the group consisting of linear, quadratic, cubic 

and butterfly. 
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35. A method as recited in any of claims 28-34 wherein the particular vertex at which 
each wavelet coefficient is defined is incident upon the midpomt of an edge belonging to two ' 
larger triangles at a next higher level of subdivision, and the wavelet coefficient is uniquely 
associated with one of the two larger triangles. 

5 36. A computer program product comprising a computer-usable medium having 

computer-readable code embodied thereon for uncompressing data representing a function 
defmed upon a two-dimensional manifold, the manifold being embedded in three-dimensional 
space, a geometric base being used as a coarse initial model of the manifold, the geometric base 
being subdivided using a triangular subdivision scheme to produce a refined triangular mesh, the 

10 computer program product comprising the following computer-readable program code for 
effecting actions in the computer 

program code for defining a tree structure in which each node of the tree structure 
represents an associated triangle of the triangular mesh, each node arranged to accommodate 
associated second generation wavelet coefficients, there being a plurality of levels in the tree, the 
15 tree levels including a top level that includes at least one top level node, wherein each tree level 
below the top level includes a plurality of nodes that are children of parent nodes in a next higher 
level, wherein children nodes associated with a particular parent node represent triangles of a 
subdivided triangle associated with the particular parent iKxle; 

program code for setting a first threshold and inputting a stream of significance bits 
20 indicative of significant nodes, wherein a significant node is defined as a node that has at least one 
of, associated wavelet coefficient that exceeds the designated threshold and descendant node that 
has an associated wavelet coefficient that exceeds the designated threshold, wherein when a 
particular significant node is determined, the input stream of significance bits further includes bits 
indicative of the significance of children nodes of the particular significant node and the input bit 
25 stream identifies the significant descendants; 

program code for inputting a stream of data bits indicative of values for second generation 
wavelet coefficients, the wavelet coefficients being defined at the vertices of the triangles in the 
refined triangular mesh, wherein each wavelet coefficient at a particular level in the triangular 
mesh is uniquely defined at a particular vertex in the triangular mesh at that level and is uniquely 
30 assigned to a single one of the triangles at a next higher level of subdivision, such that each 
triangle at the next higher level of subdivision has in the range of zero to three assigned wavelet 
coefficients; and 

program code for synthesizing the data representing the function from the wavelet 
coefficients in order to reconstruct the function defined upon the two-dimensional manifold. 
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